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ABSTRACT 

We analyse the orbital distribution of elliptical (E) and lenticular (SO) galaxies using SAURON 
integral-field stellar kinematics within about one effective (half light) radius. We construct 
the anisotropy diagram, which relates the ratio of the ordered and random motion in a galaxy 
• (V/(j), to its observed ellipticity (e), for the 48 E/SO galaxies from the SAURON survey. For 

a subsample of 24 galaxies consistent with axisymmetry, we use three-integral axisymmetric 
Q . Schwarzschild dynamical models to recover the detailed orbital distribution and we find good 

$-H ' agreement with the anisotropy derived from the {V/ cr, e) diagram. In the companion Paper IX 

I we show that the early-type galaxies can be subdivided into two classes of systems with or 

■ without a significant amount of specific stellar angular momentum. Here we show that the 

two classes have different distributions on the (y/cr, e) diagram. The slow rotators are more 
common among the most massive systems and are generally classified as E from photometry 
alone. Those in our sample tend to be fairly round (e < 0.3), but can have significant kinemat- 
ical misalignments, indicating that as a class they are moderately triaxial, and span a range of 
anisotropics {6 < 0.3). The fast rotators are generally fainter and are classified either E or SO. 
They can appear quite flattened (e < 0.7), do not show significant kinematical misalignments 
(unless barred or interacting), indicating they are nearly axisymmetric, and span an even larger 
range of anisotropics {6 < 0.5). These results are confirmed when we extend our analysis to 
18 additional E/SO galaxies observed with SAURON. The dynamical models indicate that the 
anisotropy inferred from the {V/a, e) diagram is due to a flattening of the velocity ellipsoid 
in the meridional plane (ct^ > cr^), which we quantify with the (3 anisotropy parameter. We 
find a trend of increasing (3 for intrinsically flatter galaxies. A number of the fast rotators show 
evidence for containing a flattened, kinematically distinct component, which in some cases 
counter rotates relative to the main galaxy body. These components are generally more metal 
rich than the galaxy body. All these results support the idea that fast rotators are nearly oblate 
and contain disk-like components. The role of gas must have been important for their forma- 
tion. The slow rotators are weakly triaxial. Current collisionless merger models seem unable 
to explain their detailed observed properties. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: kinematics and dynamics - galaxies: structure 
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1 INTRODUCTION 

Early-type galaxies are traditionally classified into elliptical (E) and 
lenticular (SO) galaxies. More than thirty years ago, before obser- 
vations of the stellar kinematics were available, E galaxies were 
thought to constitute a homogeneous class of simple oblate systems, 
with an isotropic velocity distribution, and in which the flatt ening 
provi ded a measure of the galaxy angular momentum (e. g. iGotd 
Lenticular galaxies were considered similar to spiral galax- 
ies, containing an isotropic bulge and a disk, but without significant 
evidence for gas and dust. 

The revolution came in the late 70's, when it became possible 
to extract the stellar kinematics of bright galaxies. To characterise 
the degree of ordered rotation in galaxies the anisotropy diagram 
was introduced, which relates the ratio between the ordered and 
random motion iV/o) in a galaxy to its observed flattening (e). The 
small observed values ofV/o, for a sample of 13 bright E galax- 
ies, was interpreted as evidence that these systems are inconsistent 
with being oblate systems, close t o isotropic and supported by rota - 
tion ( Bertola & CaDaccioli|[T97^ : llUingwortlJl 19771 : iBinnevll 19781) . 
As rotation was not needed to flatten the systems, the observations 
could also be explained by assuming that E galaxies were triaxial 
and supported by orbital anisotropy (Binney 1976, 1978). The ob- 
servations were late r extended to a sample of 11 fainter E galaxies 
toavies etal]|l983b and bulges of 9 barred aKormendvHl982h) and 
8 unbarred iKormendv & Illingworthiri982h spiral galaxies. From 
the resulting compilation of a sample of 50 Es and bulges, over two 
orders of magnitude in luminosity, it was found that, contrary to the 
brighter systems, the fainter galaxies and the spiral bulges were ro- 
tating fast and their location o n the (V/o, e) d iagram was consistent 
with oblate isotropic rotators jPavies et aljl9 83). 

It was later discovered that the fast rotating galaxies also tend 
to have di sky isophotes, while the slow ro tating ones have boxy 
isophotes ( lBendej|l988l : lBender et aljl989b . One possibility to ex- 
plain this connection between photometry and kinematics was to 
assume the fast rotating galaxies all contained disks seen at var- 
ious inclinations on top of a spheroidal stellar c omponent (e. g. 
iRix & Whitj|l99(]| : iRix, Carollo. & Freemaj[l999i) . Observations 
at different scale heights on edge-on disk galaxies however showed 
that bulges themselves are rotating fast Kormendv & lUingworthl 
[l982i). The slow rotation of the brighter Es could be equally due 
to either radial or tangential anisotropy. To address this point 
Ivan der Marell Il99ll) used two-integral Jeans dynamical modelling 
to analyse a compilation of long slit kinematics of 37 bright Es. He 
concluded that Es as a class are radially anisotropic. This appeared 
consistent with what had been observed i n the outer parts o f A'^- 
body simulations of collisionless collapse jvan AlbadjI 19821) . The 
difference between the different types of spheroidal systems could 
be explained in the context of galaxy f ormation as due to the rel- 
ative importance of gaseous dissipation jBender. Burstein. & Fabej 

Il992h . 

The whole picture was summarised by iKormendv & Bended 

j 19961) who proposed a revision of the standard photometric classi- 
fication daibble 1936) of early-type galaxies into E and SO. They 
suggested that there was a dichotomy between two classes of early- 
type galaxies: (i) on one side the bright, boxy, slow rotating, and 
radially anisotropic E galaxies, with a break between a steep outer 
surface brightness profile and a more shallow nuclear cuspy core, 
and (ii) on the other side the generally fainter, disky, fast rotat- 
ing and isotropic disky-E and SO, without clear breaks between the 
outer and nuclear profiles. 

More recently a systematic investigation of the anisotropy 



of a sample of 21 luminous, nearly round and slowly rotat- 
ing elliptical gala xies was perfor med using sph erical models by 
iKronawitter et al.l (2000) and Gerh ard et alUzGOll) . They found that 
elliptic als are only moderately radially anisotropic. iGebhardt et alj 
l l2003l) constructed more general axisymmetric models and ex- 
plored the anisotropy of a sample of 12 early-type galaxies, which 
included flattened objects. They found a range of anisotropy and no 
obvious trend (except near the nuclear supermassive black holes), 
but they also found the roundest galaxies to be generally only 
weakly anisotropic. However the models in these studies were gen- 
erally fitted to the kinematics extracted along a few long-slit posi- 
tions. 

The introduction of panoramic integral-field spectrographs on 
large telescopes, combined with the advances in the dynamical 
modelling techniques, has opened the possibility for a reanalysis of 
the orbital distributi on in early-type ga laxies. The goal of this paper 
is to use SAURON jSacon et al.ll200ll hereafter Paper I) integral- 
field stellar kinematics t o derive the (V/a, e) anisotropy diagram 
for the SAURON sample jde Zeeuw et alj|2002l hereafter Paper II) 
in a ne w way t hat ca n be rigorously interpreted with the formal- 
ism of iBinnevI ( |2005[) . The results are analysed making use of a 
new classification scheme that we introduce in a companion pa- 
per iEmsellem et al.,,2007^ hereafter Paper IX). The anisotropy de- 
rived from the {V/ a, e) diagram is interpreted and tested using gen- 
eral axisymmetric dynamical models for a subsample of the survey 
galaxies. 

The galaxy sample is described in Section 2. In Section 3 we 
describe the dynamical modelling technique and in Section 4 we 
present our modelling results. In Section 5 we present the iV/cj, e) 
diagram obtained from the SAURON data, while in Section 6 we 
compare our results with previous works. Finally in Section 7 we 
draw some conclusions. 



2 SAMPLE AND DATA 

Two different subsamples are used in this paper, (i) The main galaxy 
sample (Table [T]l is composed of the 48 E and SO galaxies of the 
SAURON survey (Paper II). The {V/o, e) diagram of this sample is 
presented in Section[5]2] (ii) For half of the survey galaxies we con- 
structed axisymmetric dynamical models to interpret the SAURON 
{V/(T,e) diagram. This mo delling subsampl e of 24 galaxies (Ta- 
ble |2j is the same used in ICappellari et al.l ( l2006l hereafter Pa- 
per IV). 

All the galaxies used have SAURON integral-field spectroscopy 
out to about one effective (half lig ht) radius (Re). T he SA URON 
data were reduced as described in IEmsellem et al. I ( l2004 here- 
after Paper III). However there are some differences between the 
kinematics presented in Paper III and the one used for this work: 
(i) To provide a tight constraint to the dy namical models we ex- 
tracted the Gauss-He rmite (G-H) moments I van der Marel & Franxl 
ll993l : lGerhitf3ll993h of the line-of-sight stellar velocity distribu- 
tion (L OS VP) up to fca-fee, u sing the penalised pixel-fitting method 
(pPXF. iCappellari & Emsellem.2004.) : (ii) To measure the mean ve- 
locity V and the velocity dispersion cr to be used in the {V/a, e) 
diagram we did not fit the higher order G-H moments. We verified 
using pPXF on semianalitic dynamical models that this generally 
provides a better approximation to the first and second velocity mo- 
ment of the LOSVD, which appear in the tensor-virial equations 
from which the {V/a, e) diagram is constructed; (iii) To reduce the 
influence of template mismatch we updated the library of templates 
which is fitted together with the kinematics by pPXF. In particu- 
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Table 1. Measured parameters for the 48 E/SO of the SAURON sample. 



Galaxy Name Type T i?max/i?o e (V^) (V/a)^ PAphot PAki„ APAki^ A^^ Rotator 









(arcsec) 






km s ^ 


km s ^ 




(deg) 


(deg) 


(deg) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


NGC 474 


SO"(s) 


-2.0 


29 


0.71 


0.11 


31 


146 


0.21 


14.0 


152.0 


4.5 


0.20 


F 


NGC 524 


SO+(rs) 


-1.2 


51 


0.61 


0.05 


64 


222 


0.29 


48.4 


40.0 


1.0 


0.28 


F 


NGC 821 


E6? 


-4.8 


39 


0.62 


0.40 


48 


182 


0.26 


32.2 


31.0 


1.5 


0.26 


F 


NGC 1023 


SBO"(rs) 


-2.7 


48 


0.56 


0.33 


62 


176 


0.35 


87.0 


89.0 


1.0 


0.38 


F 


NGC 2549 


SO"(r)sp 


-2.0 


20 


1.16 


0.49 


66 


119 


0.56 


1.2 


1.0 


0.5 


0.54 


F 


NGC 2685 


(R)SBO+pec 


-1.1 


20 


1.34 


0.62 


62 


71 


0.88 


36.6 


37.0 


1.5 


0.72 


F 


NGC 2fiQ5 


SARO'^Ts^ 


-2.1 


21 


0.96 


0.29 


87 


162 


0.54 


169.3 


175.0 


1.0 


0.56 


F 


NGC 2699 


E: 


-5.0 


14 


1.41 


0.15 


48 


1 12 


0.43 


50.2 


50.0 


3.0 


0.45 


F 


NGC 2768 


E6: 


-4.3 


71 


0.39 


0.38 


51 


209 


0.24 


94.7 


95.0 


1.5 


0.27 


F 


NGC 2974 


E4 


-4.7 


24 


1.04 


0.37 


127 


180 


0.70 


43.5 


43.0 


0.5 


0.60 


F 


NGC 3032 


SABO'^(r) 


-1.8 


17 


1.19 


0.11 


23 


87 


0.27 


94.6 


91.0 


10. 


0.42 


F 


NGC 31 56 


SO: 


-2.4 


25 


0.80 


0.47 


42 


48 


0.88 


49.4 


50.0 


3.5 


0.71 


F 


NGC 3377 


E5-6 


-4.8 


38 


0.53 


0.46 


57 


1 17 


0.49 


41.3 


46.0 


1.0 


0.47 


F 


NGC 3379 


El 


-4.8 


42 


0.67 


0.08 


28 


198 


0.14 


67.9 


72.0 


2.0 


0.14 


F 


NGC 3384 




-2.7 


27 


0.75 


0.20 


57 


131 


0.44 


53.0 


48.0 


1.5 


0.41 


F 


NGC 34 1 4 


SO pec 


-2.1 


33 


0.60 


0.21 


19 


206 


0.09 


179.9 


1.0 


9.5 


0.06 




NGC 3489 


SAR0+(rs1 


-1.3 


19 


1.05 


0.29 


52 


78 


0.67 


71.2 


73.0 


1.0 


0.60 


F 


NGC 3608 


E2 


-4.8 


41 


0.49 


0.18 




179 


0.05 


79.3 


85.0 


8.5 


0.04 


s 


NGC 4150 




-2.1 


15 


1.39 


0.30 


37 


64 


0.58 


147.0 


147.0 


4.5 


0.58 


F 


NGC 4762 


SRO" fs^ 


-2.7 


10 


2.06 


0.22^ 


40 


166 


0.24 


154.0 


150.0 


1.5 


0.25 


F 


NGC 4270 


so 


-1.9 


18 


1.09 


0.50 


46 


1 14 


0.40 


107.5 


102.0 


3.0 


0.45 


F 


NGC 4278 

1 'I V.J V... / O 


El-2 


-4.8 


32 


0.82 


0.12 


44 


228 


0.19 


16.7 


12.0 


0.5 


0.15 


F 


NGC 4374 


El 


-4.2 


71 


0.43 


0.15 


7 


282 


0.03 


128.2 


141.0 


1.5 


0.02 


s 


NGC 4382 


SO+TsWr 


-1.3 


67 


0.38 


0.19 


31 


192 


0.16 


29.4 


23.0 


1.5 


0.16 


F 


NGC 4387 




-4.8 


17 


1.16 


0.40 


34 


89 


0.39 


141.5 


151.0 


6.0 


0.41 


F 


NGC 4458 


EO-1 


-4.8 


27 


0.74 


0.12 


10 


84 


0.12 


4.5 


24.0 


17 


0.05 


s 


NGC 4459 


SO+Tr^ 


-1.4 


38 


0.71 


0.17 


66 


146 


0.45 


102.7 


100.0 


0.5 


0.44 


F 


NGC 4473 


E5 


-4.7 


27 


0.92 


0.41 


41 


188 


0.22 


93.7 


92.0 


1.0 


0.19 


F 


NGC 4477 


SBO(s):? 


-1.9 


47 


0.43 


0.24+ 


33 


158 


0.21 


64.0 


70.0 


2.5 


0.22 


F 


NGC 4486 


FO-1 +npr 


-4.3 


105 


0.29 


0.04 


7 


306 


0.02 


158.2 


1 13.2 


45. 


0.02 




NGC 4526 


SABO"(s) 


-1.9 


40 


0.66 


0.37 


103 


189 


0.54 


112.8 


111.0 


0.5 


0.47 


F 


NGC 4546 




-2.7 


22 


0.94 


0.45 


97 


161 


0.50 


75.0 


79.0 


0.5 


0.60 


F 


NGC 4550 


SBO^'isp 


-2.0 


14 


1.45 


0.61 


12 


1 16 


0.10 


178.3 


178.0 


1.0 


0.09 


S* 


NGC 4552 


EO-1 


-4.6 


32 


0.63 


0.04 


13 


257 


0.05 


125.3 


1 13.0 


5.0 


0.05 


s 


NGC 4564 




-4.8 


21 


1.02 


0.52 


76 


131 


0.58 


48.6 


49.0 


2.0 


0.59 


F 


NGC 4570 


SO sp 


-2.0 


14 


1.43 


0.60 


81 


152 


0.53 


159.3 


159.0 


0.5 


0.56 


F 


NGC 462 1 


E5 


-4.8 


46 


0.56 


0.34 


52 


207 


0.25 


163.3 


165.0 


0.5 


0.27 


F 


NGC 4660 




-4.7 




1.83 


0.44 


79 


163 


0.49 


96.8 


98.0 


0.5 


0.47 


F 


NGC 5198 


El-2: 


-4.7 


25 


0.80 


0.12 


12 


185 


0.07 


15.3 


58.0 


21. 


0.06 


s 


NGC 5308 


SO- sp 


-2.0 


10 


2.04 


0.60 


86 


192 


0.45 


58.5 


59.0 


1.0 


0.48 


F 


NGC 5813 


El-2 


-4.8 


52 


0.53 


0.15 


32 


223 


0.14 


134.5 


151.0 


2.5 


0.06 


S 


NGC 5831 


E3 


-4.8 


35 


0.67 


0.15 


11 


151 


0.08 


122.8 


101.0 


10. 


0.05 


s 


NGC 5838 


SO- 


-2.7 


23 


0.87 


0.34 


110 


216 


0.51 


41.7 


39.0 


0.5 


0.52 


F 


NGC 5845 


E; 


-4.8 


4.6 


4.45 


0.35 


81 


226 


0.36 


143.2 


141.0 


2.0 


0.36 


F 


NGC 5846 


EO-1 


-4.7 


81 


0.29 


0.07 


7 


240 


0.03 


75.2 


126.0 


5.0 


0.02 


s 


NGC 5982 


E3 


-4.8 


27 


0.94 


0.30 


19 


234 


0.08 


108.9 


114.0 


4.0 


0.09 


s 


NGC 7332 


SO pec sp 


-2.0 


11 


1.91 


0.42 


38 


116 


0.32 


159.8 


152.0 


1.5 


0.39 


F 


NGC 7457 


SO-(rs)? 


-2.6 


65 


0.39 


0.44 


38 


62 


0.62 


125.5 


124.0 


4.0 


0.57 


F 



Notes: (1) NGC number. (2) Morphological type from lde Vaucouleurs et aljjl99 1. hereafter RC3). (3) Numerical morphological T-type (LEDA. E: T < —3.0, 
SO: -3.0 < T < -0.5). (4) Effective (half-light) radius measured in the 7-band from HSTAVFPC2 + MDM images as described in Paper IV. Comparison 
with the RC3 values, for the 46 galaxies in common, shows an rms scatter of 20%. (5) Ratio between the maximum radius iJmax sampled by the kinematical 
observations and R^. We defined ifmax = \/ S/ir, where S is the area on the sky sampled by the SAURON observations. (6) Luminosity-weighted average 
ellipticity. This was computed from the ellipse of inertia of the galaxy surface brightness inside an isophote enclosing an area A = wRc'^, or within the largest 
isophote fully contained within the SAURON field, whichever is smaller. (7) Luminosity-weighted squared velocity within an ellipse of area A, ellipticity e, 
and PA given in column [10], or within the largest similar ellipse fully contained within the SAURON field, whichever is smaller. (8) Luminosity-weighted 
squared velocity dispersion inside the same ellipse as in column [7]. (9) Luminosity- weighted V/cr ratio within liJc. This is the ratio of columns [7] and [8]. 
See equation 1101 for a definition of these quantities. (10) Large scale global luminosity-weighted PA of the photometric major axis. (11) PA of the global 
kinematic major axis within the SAURON field (direction where IV] is maximum, see text for details). (12) Error in the kinematic major axis of column [11]. 
(13) Afl (see Paper IX) measured within the same ellipse as in column [7] (about IRc)- (14) Galaxy classification from Paper IX: F=fast-rotator (Ajij^ > 0.1), 
S=slow-rotator (Aj:j < 0.1). 

* This galaxy is a special slow-rotators. It appears axisymmetric and disk-like as a fast-rotator, but contains two counterrotating disks (Section [4.5t . 
t These two galaxies show nearly face-on bars. The ellipticity is the one of the outer disk, which also defines the photometric PA. 
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Table 2. Anisotropy parameters for the 24 modelled galaxies. 



Galaxy Name 



(deg) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




1 Q 


n r\A 
U.Uo 


n 1 7 
U. 1 / 


n r\A 
-U.U4 


n 1 Q 

u. ly 


n 1 o 

u. ly 




on 


n 1 ^ 
u. io 


n 7 1 
U.Zl 


n n/1 

U.U4 


n 7n 
u.zu 


n in 
U. JU 


Mr;r^ TQ7/I 

IN Kjy^ Ay 1 


^7 


n on 


U. 1 J 


n in 

-U. jU 


n 7/1 

U.Z4 


n 7n 
u.zu 




Do 


n 1 7 


n 70 

V.Jy 


n 1 Q 
u. ly 


nil 

U.J J 


n 1 
u. iy 


Mnr^ 1177 


on 


n n7 


n 7Q 
U.Zo 


n ns 
U.Uo 


n 7^ 

U.Zj 


n 71 

U.Zj 


Mnr^ 1170 


on 


n 1 1 

U. i 1 


n n*^ 

U.UO 


n nA 

U.UO 


n m 
u.u J 


n r\A 

U.U4 




on 


n 1 1 


n n*^ 

U.UO 


n 1 7 
-U. iz 


n 1 1 

U. 1 i 


n 1 7 

U. i / 




on 


n r\A 

U.U4 


n 1 n 

U. iU 


n n*^ 

-U.UO 


n 1 1 

U. 1 J 


n 1 ^ 

U. i J 




JZ 


n ni 

-U.Ui 


n T7 

U. jZ 


n 1 7 
-u. iz 


n iA 
U. jO 


n 77 
u.zz 


Mnr^ /n7Q 

IN OL^ 4Z / o 


on 


n m 

-U.UZ 


n 1 1 

U. 1 1 


n 1 7 
-U. 1 / 


n 1 fi 

U. 1 o 


n 

U.UO 


Mr;r^ /1 17/1 

IN LjL^ 4-5 / 4 


on 
yu 


nil 

U. i 1 


n 1 n 

U. iU 


n n^ 

U.UJ 


n ns 

U.Uo 


n 1 7 

U. iZ 




on 
yu 


n 7A 
-u.zo 


n ni 

-U.Ui 


n 71 


n no 
u.uy 


n nc 

U.Uo 




Al 


n 1 n 


n ns 

U.UJ 


n 1 1 

U. 1 1 


A nn 
u.uu 


n n'? 

U.UZ 


NGC 4473 


73 


-0.21 


0.18 


-0.50 


0.34 


0.37 


NGC 4486 


90 


0.24 


0.11 


0.22 


0.00 


0.03 


NGC 4526 


79 


0.11 


0.11 


0.09 


0.06 


0.09 


NGC 4550 


84 


-0.37 


0.43 


-0.87 


0.60 


0.56 


NGC 4552 


90 


-0.06 


0.01 


-0.03 


0.02 


0.03 


NGC 4621 


90 


-0.04 


0.11 


-0.17 


0.18 


0.24 


NGC 4660 


70 


0.02 


0.27 


-0.11 


0.30 


0.30 


NGC 5813 


90 


0.17 


0.18 


0.21 


0.08 


0.10 


NGC 5845 


90 


0.24 


0.23 


0.18 


0.15 


0.19 


NGC 5846 


90 


0.17 


0.09 


0.17 


0.01 


0.06 


NGC 7457 


64 


0.03 


0.38 


0.04 


0.37 


0.31 



Notes: (1) NGC number (2) Inclination from paper IV. (3) anisotropy 
parameter f3r measured in spherical coordinates from the solution of the 
dynamical models as defined in equation (8). (4)-(6) anisotropy param- 
eters 13, 7 and 5 determined in cylindrical coordinates from the solution 
of the dynamical models. The parameters are defined in equations I4I5I 
and equation (2) respectively. (7) Anisotropy parameter S as measured 
using the {V/cr, e) diagram. 



lar we construc ted the optimal template using the 985 stars of the 
MILES library ISanchez-Blazqu ez et al. from which ~ 15 

stars are selected by the program to provide a detailed fit to each 
galaxy spectrum. The use of this new library can reduce the rms 
scatter in the residuals of the pPXF fit by up to a factor ~ 3, when 
the spectra have negligible Poissonian noise. One can expect a re- 
duction of the systematic errors in the G-H moments by up to the 
same factor (see Fig. B3 of Paper III). In practice our results are 
very similar to the ones presented in Paper III, but for some of the 
most ma ssive galaxies the h 4 values are here significantly lower 
(see also lShapiro et al.ll2006l) . 

Together with the HSTAVFPC2 photometry in the J-band, 
wide-field ground-based MDM photometry (Falcon-Barroso et al. 
in preparation) is also available for all the modelled galaxies. We 
used this set of photometric data to parameterise the stellar density 
distribution in our m odels according to the multi-Gaussian e xpan- 
sion (MGE) method ( Emse llem et"5]| 1994 Ic"appellarill2002 ). The 
parameters of the MGE models are given in Paper IV. The MDM 
photometry was also used to determine the Rc of the sample galax- 
ies. 

An additional set of 18 E/SO galaxies has been independently 
observed with a similar SAURON setup in the course of various 
other projects. These objects will be treated as "specials", most 
having features which motivated a specific observation, and will 
be only mentioned in Section lS!4l to strengthen the results obtained 
from the main survey. As in Paper II the galaxies are classified as 



either E or SO fro m the LEDA morphological type (T < —0.5; 
IPaturel et all2003h . The kinematics of these galaxies were not pre- 
sented in Paper III, and will be presented elsewhere. 



3 THREE-INTEGRAL DYNAMICAL MODELLING 

The models we study in this paper were presented in Paper IV. 
The stationary and axisymmetr ic dynamical m o dellin g technique 
that we use is based on the ISchwarzschildl ( 1 19791) numerical 
orbit-superpositio n method, which was extended to fit kinemat- 
ical observables jRichstone & Tremaind 1 19881: iRix et alj 1 19971 : 
Ivan der Marel et alj|l998h . The implementation of the method that 
we adopt in this paper was optimised for use with integral-field data 
and is described in Paper IV. Similar axisymmetric implementations 
were develo ped by other groups (Gebhardt et al. 2003; Valluri et alJ 
l2004lThomas et al]|20oA . 

For a given stationary gravitational potential, the stellar dy- 
namics of a galaxy is uniquely defined by the orbital distribution 
function (DF), which describes the velocities of the stars at every 
position in the galaxy. As the stellar orbits in a stationary poten- 
tial conserve at most three isolating integrals of motion, the DF 
can be written as a function of the three integrals, or any other 
equivalent parameterization of them. From dimensionality argu- 
ments this implies that the three-dimensional DF cannot be re- 
covered without at least the knowledge of the LOSVD at every 
spatial position {x',y') on the galaxy image on the sky, which 
also constitutes a three-dimensional quantity. For an axisymmet- 
ric edge-on galaxy, with a given potential, this knowledge of the 
LOSVD seems likely sufficient for a unique recovery of the DF, 
which may not be positive everywhere if the assumed potential 
is wrong (see Section 3 of IValluri et alj |2004|) . An example of 
the need for two-dimensional kinematics to constrain the DF is 
given in Cappellari & McDermid (2005). In Krajnovic et al. (2005) 
and lvan de Ven. de Zeeuw. & van den BoschI ( 120071) we verified the 
ability of our implementation of the Schwarzschild method to re- 
cover the DF and the internal velocity moments for both two and 
three-integral realistic galaxy model, even from incomplete radial 
coverage. 



An implicit assumption of the dynamical models is that the 
luminosity density, as can be obtained by deprojecting the galaxy 
surface brightness, provides a good description of the shape of the 
total density. This is also a key assumption of the {V/a, e) dia- 
gram. It implies that either dark matter provides a small contribu- 
tion to the total matter in the regions we study (inside IRa), or that 
its shape is similar to that o f the luminous matter. Evidenc e from 
dynam ical modelling (e.g. iGerhard et alj [2 001 ; Cappell ari et alj 



2006h and g ravitational lensing (e.g. Tre u & KoopmansI 2004 



Rusin & Kochanek 200^ ; lKoopmans al. 200® suggest that both as 



sumpti ons are reasonably well Justified in real galaxies. 

In iKrainovic et al. I ( l2005l) and in Paper IV we showed that 
there is evidence for the inclination to be possibly degenerate, even 
with the knowledge of the LOSVD at all spatial positions. As dis- 
cussed in Paper IV, we adopted for our models the assumption- 
dependent inclination derived by fitting two-integral Jeans models. 
This inclination appears to provide values in agreement with the 
geometry of dust/gas disks when they are present. However for 14 
of the 24 modelled galaxies the inclination is already constrained 
by arguments independent from the dynamics. In fact nine galax- 
ies (NGC 821, NGC 3156, NGC 3377, NGC 4473, NGC 4550, 
NGC 5845, NGC 4621, NGC 4660, NGC 5845), show significantly 
disky isophotes and require Gaussians flatter than q' < 0.4 in their 



© 2007 RAS, MNRAS OOP. [711271 



The SAURON project -X 5 



MGE models (see Paper IV). Under the assumption of axisym- 
met ry, this implies t hey cannot be too far from edge-on (i > 70; 
see ICapDellari[l2002l § 2.2.2). The inclination of five of the re- 
maining galaxies can be derived from the geometry of a gas or 
dust disk, assuming it is in equilibrium in the equatorial plane of 
an oblate galaxy (NGC 524, NGC 2974, NGC 4150, NGC 4459, 
NGC 4526), and it agrees with the dynamically determined inclina- 
tion (Paper IV). 



4 MODELLING RESULTS 
4.1 Internal velocity moments 

Here we present the orbital distribution of the 24 early-type galax- 
ies of our modelling subsample (Section derived from our 
Schwarzschild models, at the inclination determined in Paper IV. 
The models provide a full description of the DF in terms of a three- 
dimensional set of weights, which could in principle be transformed 
into a distribution of stellar mass as a function of the three integrals 
of motion (E, Lz, and a nonclassical third integral I3). In this pa- 
per we will mainly focus on a more concise description of the main 
characteristics of the internal kinematics which can be obtained by 
computing the first and second moments of the velocity distribu- 
tion. Given the assumed axial symmetry, it is sufficient to study the 
moments in the galaxy meridional plane {R, z). 

The results of this section are based on a sample of galaxies 
which was selected to be consistent with the axisymmetric assump- 
tion. The sample contains both fast-rotating and nearly non-rotating 
galaxies, but it does not include strongly triaxial ones. This has to 
be considered when drawing statistical conclusions from the mod- 
els. The results of this section are important to interpret and test 
the more general results we will obtain using the {V/a, e) diagram 
in Section |5] for the whole SAURON sample which also includes 
significantly triaxial galaxies. 

A large number of studies, in the past twenty years, has shown 
that elliptical galaxies display a dichotomy between the most mas- 
sive galaxies, which tend to rotate slowly, are metal-rich, have a 
flat central luminosity profile and show evidence of triaxiality, and 
the less massive objects, which rotate faster, are metal-poor, have a 
cuspy luminosity profile, and may all b e axisymmetric and contain 
disks j Pavies et al. 1983; Bender 1988'; lKormendv & BendeJig?^ : 
iFaber etal. 1997 ; Lauer et al. 2005) . 

In the companion Paper IX we show that a clean separation of 
the two classes of early-type galaxies is obtained by introducing a 
new quantitative classification parameter Xh which is related to the 
speci fi c angular momentum of galaxies (e.g. lBertola & Capacciolil 
1 19781 ; iBender & Nietd [l990i) and which can be measured from 
integral-field stellar kinematics as 

^ ^ {R\v\) ^ 

where Fn is the flux contained inside the n-th Voronoi bin and Vn 
and a„ the corresponding measured mean velocity and velocity dis- 
persion (see Paper IX). The numerical values of Xr for our sample, 
measured with IRc, are given in Table[T] In this work we adopt that 
classification to define the galaxies in our sample with significant 
angular momentum per unit mass Xr > 0.1, which we call "fast- 
rotators" and the ones with negligible amount of specific angular 
momentum < 0.1, which we define as "slow-rotators" (Ta- 
ble [TJ. Our kinematic classification has a similar spirit to the one 
proposed by Kormendv & Bender. (,1996i) . but contrary to the latter. 



which was based on the galaxy isophotal shape, our classification is 
quite robust against projection effects. 

In Fig.[T]we plot the internal moments in the meridional plane 
of our subsample of 24 modelled galaxies. The ellipses represent 
the cross section of the velocity ellipsoids at various positions in 
the galaxy body, while the colours visualise the ratio {v^)/<7 be- 
tween the mean velocity around the symmetry axis and the lo- 
cal mean velocity dispersion a'^ = (cr^ + + a'^)/3, where 
""0 ~ (^0) ~ {'"4)'^ ■ Here (r,9,(j)) are the standard spherical co- 
ordinates ((/!> is the azimuthal angle around the galaxy symmetry 
axis). The quantity {v^}/a is a local analogue, inside the galaxy 
body, of the global quantity V/a which measures the importance of 
rotation with respect to the random motions. In the plot the galaxies 
were sorted according to the parameter. The first nine galax- 
ies are classified as slow-rotators, while the remaining fifteen are 
fast-rotators. 

The slow rotators in this modelling subsample are all intrin- 
sically quite round, with an axial ratio of the density q > 0.8. This 
result depends on the correctness of the assumed inclination, which 
in all cases is best fitted as edge-on (i — 90°). This inclination can- 
not be correct in all cases, however these objects are all very round 
in projection, which statistically implies they are intrinsically not 
far from spherical. The flattest slow-rotator of this subsample has 
an apparent axial ratio of q' ~ 0.8 and assuming it has the mean 
inclination of a randomly oriented sample (i — 60°), it would still 
have an intrinsic axial ratio q ~ 0.7. As we verified numerically 
with our models, for nearly round models the results for the or- 
bital distribution change little with inclination, as can be understood 
from the obvious fact that a spherical model appears identical from 
any inclination. In particular we constructed models for NGC 3379, 
NGC 4486 and NGC 4552 at an inclination of t = 45° and found 
virtually unchanged anisotropy results. The slow rotators of this 
subsample are generally characterized by a nearly round velocity 
ellipsoid in the meridional plane and do not show any significant 
rotation. A clear exception is NGC 4550, which shows internal ro- 
tation and a strongly flattened velocity ellipsoid. 

Our anisotropy results could be significantly in error if the 
slow rotators of our modelling subsample were intrinsically flat- 
tened and nearly face-on objects (i ^ 0). However the small ap- 
parent ellipticity of the slow-rotators is also generally true for the 
whole SAURON sample. That sample includes some flatter and 
clearly triaxial slow-rotators, however the smallest apparent axial 
ratio of any slow-rotator is still q' « 0.7. Statistically this implies 
the ratio of the shortest and longest axis of the intrinsic density is 
c/a > 0.7. In Paper IX we show that the slow rotators are truly dif- 
ferent from the fast-rotators and do not simply appear different due 
to projection. All this makes it very unlikely that any of the slow- 
rotators we modelled is a flat system seen nearly face-on and im- 
plies our anisotropy results for the subsample are reliable. 

The fast rotators appear to span a large range of intrinsic flat- 
tening with q ~ 0.3 — 0.9. For the flat objects, the maximum 
contribution to rotation generally occurs on the galaxy equato- 
rial plane, as expected, but an important exception is NGC 4473, 
which has a minimum of the velocity contribution on the major 
axis (as NGC 4550). The orbital structure of the two special cases 
NGC 4473 and NGC 4550 will be explained in Section |43] The 
fast-rotators are generally significantly anisotropic, with large vari- 
ations in the ratio ar/ae. As a first approximation the velocity el- 
lipsoid tends to be aligned in polar coordinates and appears to be 
flattened in the same direction as the density distribution, in the 
sense that ar < 00 along the galaxy rotation axis, while cr,. > ae 
on the equatorial plane. Using cylindrical coordinates in the merid- 
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Figure 1. Internal moments of the modelling subsample. The colours visualise the ratio {V) /a of the first and second internal moment of the velocity, in the 
positive quadrant of the meridional plane (i?, z), for the 24 galaxies of the modelling subsample. The contours of the intrinsic density are overlaid. The ellipses 
are cross sections of the velocity ellipsoid with the (tifl, fz) plane at different locations (one of the principal axes of the ellipsoid is orthogonal to the plane). 
The objects are sorted, from left to right and from top to bottom, according to the parameter (Table[T) which is related to the specific angular momentum 
per unit mass. The first nine galaxies are classified as slow-rotators (A^^ < 0.1) and the remaining fifteen as fast-rotators > 0.1). All panels have the 

same size, which roughly corresponds to the area fuUy covered by the observed SAURON kinematics. 



ional plane [R, z), one may broadly describe the observed shape of 
the velocity ellipsoid as being generally flattened along the z direc- 
tion. There seems to be a tendency for the anisotropy to be stronger 
near the equatorial plane, especially in disk-dominated galaxies like 
NGC 3156, NGC 4150 and NGC 7457. 



4.2 Global velocity dispersion tensor 

In the previous section we described the variations of the velocity 
ellipsoid as a function of position inside the galaxies. In this section 
we consider global integrated quantities for each galaxy. A classic 
way to quantify the global anisotropy in galaxi es is by using the 
anisotropy parameter ( iBirmev & Tremain3 19871 § 4.3) 



n ' 



(2) 



where z coincides with the symmetry axis of an axisymmetric 
galaxy, x is any fixed direction orthogonal to it and 



pol d^x, = Mn (jj., 



(3) 



with (Tfe the velocity dispersion along the direction fc at a given 
location inside the galaxy. The summation defines how we com- 
puted this quantity from our Schwarzschild models. Mn is the mass 
contained in each of the polar bins in the meridional plane of 
the model, and ak,n is the corresponding mean velocity dispersion 
along the direction fc. At any location in the galaxy, the velocity el- 
lipsoid is defined by having the principal axes along the directions 



which diagonalize the tensor a. We define two additional anisotropy 
parameters: 



(4) 



describes the global shape of the velocity dispersion tensor in the 
(vn, Vz) plane. /3 = if the intersection of the velocity ellipsoid in 
the {vr, Uz) plane (as shown in Fig.[Tll is everywhere a circle. This 
is the case e. g. if the DF depends only on the two classical integrals 
of motion / — f{E,Lz), where E is the energy and Lz is the 
angular momentum parallel to the z-axis. If the shape of the velocity 
ellipsoid is constant inside the galaxy body then fJ = 1 — (ctz/itb)^. 
The second parameter 



7 ' 



1 - 



n 



4"P 



(5) 



describes the global shape of the velocity dispersion tensor in a 
plane orthogonal to Uz . 7 = when the intersection of the velocity 
ellipsoid with a plane orthogonal to Vz axis is a circle everywhere. 
I3 — "/ — S — for an isotropic system (spherical velocity ellipsoid 
everywhere). Integrating over the azimuthal angle one finds 

„ TIrr + 11^0 

i-i-xx = , (6) 

so that the three anisotropy parameters are related by 

2-7 

In the case 7 = the anisotropy can be measured directly in the 
meridional plane and the simple relation /3 ~ 5 applies. 



(7) 
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We evaluated these anisotropy parameters from the solution of 
the Schwarzschild models, restricting the volume integral of equa- 
tion 1I3J only within the radius (R < 25") fully constrained by the 
kinematics which is shown in Fig. [T] The results are shown in the 
top two panels of Fig. |2] In general we find that 5 ^ (3 (similar t o 
the non-classic two-integral models of lOehnen & Gerharjl 19931) . 
while there is no obvious trend of 7 with increasing anisotropy 5. 
This means that the observed anisotropy is mainly due to a flatten- 
ing of the velocity dispersion tensor in the z direction. However 
most of galaxies have 7^0, indicating mild radial anisotropy. Just 
a few galaxies have instead tangential anisotropy (7 < 0), and in 
particular NGC 4473 and NGC 4550 stand out (see Sectionl43]l. 

The near isotropy of the roundest galaxies is not surprising in 
this plot, in fact the anisotropy parameters are computed in cylin- 
drical coordinates [R, 2), as best suited for oblate bodies. By con- 
struction all the anisotropy parameters are zero by symmetry in the 
spherical limit. This is a known feature of the {V/a, e) diagram as 
well, where spherical objects necessarily have zero anisotropy 5. To 
measure the anisotropy of nearly spherical galaxies one can define 
a complementary anisotropy parameter, in spherical coordinates: 



with 



(8) 



(9) 



and (r, 9, (f>) the standard spherical coordinates. In the spherical 
limit, assuming the galaxy is non rotating, Ilea = = Tlrr 

by symmetry. The parameter is /3r = for an isotropic galaxy and 
is positive (negative) when the luminosity-weighted average disper- 
sion along the radial direction is larger (smaller) than the average 
dispersion along any direction orthogonal to it. The plot of jSr ver- 
sus 5 is shown in the third panel of Fig. [2] It confirms that the small 
anisotropy of the slow-rotators in the top two panels is not an arti- 
fact of the adopted coordinate system. This can also be seen quali- 
tatively in Fig. [T] which shows that the intersection of the velocity 
ellipsoids of the slow-rotators with the {vr,Vz) plane is nearly a 
circle (ere ~ (Jr). None of the roundish slow-rotators appear to be 
dominated by radial anisotropy in the regions that we observe and 
the inferred values of | ^ 0.2 imply that the average shape of the 
velocity ellipsoid for the roundish galaxies does not differ by more 
than « 10% from a spherical shape. 

The results in this section quantify what we discussed qual- 
itatively in Section [43] namely the fact that the slow-rotators in 
the modelling subsample of 24 galaxies tend to be nearly isotropic, 
while the fast-rotators span a large range of anisotropies. As a con- 
sistency test, in Appendix |A] we use two-integral Jeans models to 
show that the slow-rotators are well reproduced by isotropic mod- 
els in their central regions, while the flat fast-rotators are not. The 
same considerations on the bias of the selection of this modelling 
subsample against strongly triaxial galaxies, that we discussed in 
Section l4~n also apply to the results of this section. 

The anisotropy we measured was derived by fitting kine- 
matical data which have limited spatial coverage, so the values 
are not nec essarily representative of the galaxy as a whole. In 
iKrainoviTe t al. (2005) we showed that the orbital distribution in- 
ferred from our orbital-superposition models, in the case of integral- 
field data with limited spatial coverage, is accurate in the region of 
the galaxy for which we have kinematical data, which in our case 
generally corresponds to l_Rc. For most of the slow-rotators how- 
ever, we generally do not reach l_Re with our SAURON data. This 
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Figure 2. Global anisotropy from the dynamical models. Comparison be- 
tween the anisotropy 5 determined from the 24 Schwarzschild models of 
Fig. [T] and the anisotropy pai'ameter /3 {Top Panel), 7 {Second Panel) and 
f3r {Third Panel). Bottom Panel: Comparison between 5 from the models 
and <5jyy^j from the {V/a, e) diagram of Fig. [3] after correction for the 
inclination. The red and blue labels refer to the slow and fast rotators re- 
spectively. 
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has to be kept in mind wiien comparing our results to the results 
obtained e.g. from numerical simulations. 



4.3 Comparison with previous modelling results 

The modelling results discussed in this paper are the first based on 
integral-field stellar kinematics, which is critical for the recovery 
of the orbital distribution of flattened galaxies. No previous results 
are available using a comparable technique. In the spherical limit 
however observations along a single slit position are in principle 
sufficient to uniquely recover the DF (Dejonghe & Merritt 1992). 
The anisotropy we derive for the roundest galaxies can be compared 
with previous results obtained using spherical models. 

The largest set of spherical model s of galaxies in common 
with our sample is the one presented in ' Kronawitter et 
and analysed in Gerhard et al. (2001). The five galaxies in com- 
mon with our modelling subsample of 24 galaxies are NGC 3379, 
NGC 4278, NGC 4374, NGC 4486 and NGC 5846. The galaxy with 
the largest apparent ellipticity is NGC 4374, which has e = 0. 15, so 
our models, which take the small nonzero ellipticity into account, 
do not differ too much from spherical models, and we can com- 
pare the results obtained with different kinematical data and differ- 
ent methods. For four of the galaxies we find that the small radial 
anisotropy j3r ~ 0.1 — 0.2 that we measure does not differ by more 
than A|3t. < . 1 from the /Smcan (kindly provided by O. Gerhard) of 
lGerhardetal]j200lh . which is the unweighted mean of their radial 
anisotropy within 17?e (their Fig. 5). This error corresponds to a dif- 
ference of just 5% on the mean flattening of the velocity ellipsoid 
in these galaxies. This value likely represents the realistic system- 
atic error one can expect from this type of measurements with cur- 
rent data and techniques. For one galaxy, NGC 4278, they measure 
a strong anisotropy /Smcan ~ 0.3, while we derive near isotropy 
I3r ~ 0. The diffe rence could be explained by the fact that the data 
of I Kronawitter et al. (2000) for this galaxy are not very extend in 
radius and possess the lowest quality of the considered 

test sample. 

Another systematic study of the anisotropy of a sample of 12 
early-type galaxies, using an axisymmetric Schwarzschild method 
si milar to the one adopt ed for the present paper was performed 
bv lGebhardtetal.l ( l2003h . We can compare the results for the five 
galaxies in common that we both modelled as edge-on NGC 821, 
NGC 3377, NGC 3379, NGC 3608 and NGC 5845. For this we 
computed from our models an approximation to their mean ratio 
((7,./(7t)^^y^ (their Table 1 and 3). Contrary to our anisotropy pa- 
rameter /3r, their ratio is not directly related to the shape of the ve- 
locity ellipsoid, as it includes both random and ordered azimuthal 
streaming (see their Section 4.7). In particular according to their 
definition, unless the galaxy is spherical, an isotropic model will 
have (ar/(Jt) < 1, while a ratio of unity indicates radial anisotropy 
(Pr > 0). Trying to closely reproduce their measurements on our 
models we obtain a ratio {o'r/o't)fi ~ 1.0 for all five galaxies, 
with maximum differences of 10% Ijetween ours and their results. 
This consistency test is not as useful as the previous one, as it does 
not allow us to test our derivation of the shape of the velocity el- 
lipsoid in flattened galaxies. However for the two nearly round ob- 
jects NGC 3379 and NGC 3608 both sets of models show that these 
galaxies are very nearly isotropic. 

In summary, the comparison between our anisotropy determi- 
nations and those of previous authors, on two limited samples of 5 
models each, indicates a general agreement and suggests that our 
systematics errors are on the order of 5% on the shape of the veloc- 
ity ellipsoid. 



4.4 Comparison with the ( V/a, e) diagram 

In the previous sections we studied the anisotropy in galaxies by 
direct inversion of the observed kinematics and photometry, using 
axisymmetric dynamical models. In this section we revisit the pre- 
vious findings with the classic {V/o, e) diagram ( Binnevlll978l) . as 
was done for the last thirty years. This will provide a link between 
our results and the classic ones, and in addition it will allow us to 
test the robustness of the derived values. 

In Fig. [3] we visualise the location of the slow and fast rota- 
tors samples in the {V/a, e) diagram. Traditionally the observed 
V/ a quantity is computed from the central velocity disp ersion and 
the maximum rotational velocity (but see iBacorj [ 19850 . Recently 
iBirmevI ( |2005| ) updated the formalism to compute the quantity in a 
more robust way, from integral-field data. Here the availability of 
SAURON integral-field kinematics allow us to apply this new for- 
malism for the first time to real galaxies. We use the updated for- 
mulae and define 



(10) 



as a luminosity-weighted quantity, which we estimate from the 
binned SAURON kinematics. Here is the flux contained inside 
the n-th Voronoi bin and V„ and (t„ the corresponding measured 
mean velocity and velocity dispersion. Similarly we define the el- 
lipticity as 



(1 



V F 



(11) 



where the (x, y) coordinates are centred on the galaxy nucleus and 
the X axis is aligned with the galaxy photometric major axis. We es- 
timate e from the individual SAURON pixels, inside a given galaxy 
isophote, within the same region used for the computation of V/a. 
This way of measuring e by diagonalizing the inertia tensor of 
the galaxy surface b rightness is the stan dard technique among A'^- 
body modellers (e.g. lCarter & Metcalfell98Q) an d for the automatic 
determ ination of e in large galaxy surveys (e.g. iBertin & AmoutsI 
Il996h . When the ellipticity varies little with radius the e mea- 
sured with this tech nique is very similar to the intensity-weighted e 
dRvden et al.lll999h from the standard photometry profiles (see also 
Fig.lI4b. In general however the e determinations from the moments 
are biased towards larger radii than e determinations from the pro- 
files. 

For a rigorous application of the iV/a, e) diagram, the sum- 
mation JlOb should extend to infinite radii and the e should be con- 
stant with radius. In Appendix iBl we show that reliable results for 
the {V/a, e) values (at least in the limited isotropic case considered) 
can still be obtained when the summation is spatially limited and the 
ellipticity is slowly varying. To make the measurement comparable 
between different galaxies, we limit the summation to \Rc, or to 
the largest observed radius, whichever is smaller. We also show that 
better results are obtained when the summation is performed inside 
ellipses instead of circles. For this reason the quoted values of V/o 
are measured inside ellipses of area A = -kRb^ and semi-major axis 
a — Rc/\/l — e, where the mean ellipticity e is measured inside 
an isophote of the same area A. When the isophote/ellipse are not 
fully contained within the observed SAURON field of view, we limit 
the summation to the largest isophote/ellipse contained within the 
field. The {V/a, e) values are given in Table[T] 

The {V/ a, e) diagram for the 24 galaxies of the modelling sub- 
sample is shown in Fig. [3] The solid curves in the {V/a, e) dia- 
gram represent the location of oblate edge-on models with differ- 
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Figure 3. {V/a, e) diagram for our modelling subsample of 24 galaxies. 
The red and blue labels refer to the NGC number of the slow and fast rota- 
tors respectively, and show the measured values of the luminosity weighted 
ellipticity e and V/cr. The soHd lines, starting from each object show the 
effect of correcting the observed values of each galaxy to an edge-on view 
(diamonds). For the fast rotators we adopted the inclination derived from 
the models of Paper IV. For the slow-rotators we show the correction for 
inclination assuming the average inclination of a randomly oriented sam- 
ple i = 60°. Even under this assumption the corrections are small for the 
slow-rotators. The grid of solid curves shows the location on this diagram 
of edge-on oblate galaxies with different anisotropy S = 0, 0.1, ... , 0.6 
(edge-on isotropic models <5 = are shown with the thick green line). 



ent anisotropy 5. Real galaxies are generally not seen edge-on, so 
a correction for the effect of inclination is in principle needed to 
correctly interpret the diagram. This usually cannot be done, as the 
inclination of early-type galaxies is not easy to determine. In our 
case however, we have an estimate for the inclination of our galax- 
ies (Section [3), so we can correct the position of our objects on 
the diagram to the edge-on case. Although the inclination we use 
is assumption-dependent, we showed in Paper IV that it generally 
provides a big improvement from a complete lack of knowledge of 
the inclination. 

For a galaxy observed at an inclination i (edge-on correspond- 
ing to j = 90°), with an ellipticity e and a given ratio (V^/o')obs. the 
values corrected to an edge-on view are ( jsirmev & Tremain3l 19871 , 
§4.3): 



V\ VI - Sec 
obs sin i 



Einti 



= l-v/l + e(e-2)/f 



(12) 



(13) 



Starting from a given set of edge-on par ameters the an isotropy un- 
der the oblate assumption is obtained as ( lBinnevll2005h 



(V/a) 



[1 - a(V/a)2] n(e) 

where 

e = \/l - (1 - Eintr)^, 



(14) 



(15) 



0.5 (^arcsin e — e\/l — e-^) 



(1 



(16) 



and Q is a dimensionless number, which does not depend on the 
galaxy amount of rotation, but only on how the stellar density p 
and streaming velocity are distributed in the {R, z) plane. In 
Appendix IbI we use Jeans models to show that a value a ~ 0.15 
places isotropic models, with a realistic surface brightness distribu- 
tion, close to the isotropic line on the {V/o, e) diagram. Although a 
may in principle vary for different galaxies, tests suggest variations 
to be small and we will use this value for all the (V/o, e) diagrams 
in this paper. As the inclination correction requires the knowledge 
of both the known inclination i and the unknown anisotropy 5, 
one has to proceed in an iterative manner. However this process 
converges very rapidly, as the correction depends very little on 5. 
The effect of correcting the observed values of our galaxies on the 
[V/o, e) diagram is shown in Fig.|3] 

The comparison between the anisotropy derived from the 
inclination-corrected {y/o,e) diagram and from the dynamical 
models is shown in the bottom panel of Fig. [2] The generally good 
agreement confirms the reliability of the two approaches and the 
significant anisotropy of most of the fast rotators in the modelling 
subsample. 

As poin ted out by iKormendvl il982al) . in the formalism of 
iBinnevI ( 1 197 8) the isotropic line in the {V/a, e) diagram is approx- 
imated to good accuracy by the expression {V/aY ~ e/(l ~ £)■ 
Similarly, in the revised formalism for integral-field kinematics of 
Binnev (2005), which we adopt in this paper, the minimax rational 
approximation of order (1,1) to the isotropy line (equation II14I with 
(5 = and a = 0) is given by 



(17) 



which has a maximum relative error of just 0.3% over the whole 
interval e = [0.01, 0.99]. For the adopted value of a = 0.15 (Ap- 
pendix|Bll the corresponding approximation to the isotropic line be- 



0.831 



0.896 e' 



(18) 



which has a maximum fractional error of just 0.2% in the same 
interval. 



4.5 Understanding the anisotropy 



We saw in Section 1431 that the anisotropy of early-type galaxies 
tends to be due to a flattening of the velocity ellipsoid along the 
z direction, parallel to the symmetry axis. However two galaxies, 
NGC 4473 and NGC 4550, stand out for being dominated by tan- 
gential dispersion. Understanding what makes some object special 
is key to understand the normal galaxies. 

In Fig.|4]we show the data versus model comparison for the 
galaxy NGC 4550, while the same comparison for NGC 4473 was 
shown in Fig. 2 of ICappellari & McDermid j2005l) . Along the ma- 
jor axis both galaxies show a peculiar decrease in the V field and 
a corresponding increase in t he a field. As already poi nted out 
by Rubin, Graham, & Kenne^j ( Il992l) and lRix et al] (l992h , for the 
case of NGC 4550, a natural explanation for this behaviour is the 
presence of two counterrotating stellar components, as indicated by 
the observed double-peaked line profiles. In Fig.[5]this explanation 
is confirmed from the solution of the Schwarzschild models, which 
shows two major kinematically distinct stellar components, rotating 
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Figure 4. Data-model comparison. Each column from le ft to right shows the kinematical moments: mean velocity V , velocity dispersion a, and higher 
Gauss-Hermite moments hs-hg. Top row: Voronoi binned iCappellari & Copi3 l2003) and bi-symmetrised SAURON observations. Bottom row: Three-integral 
Schwarzschild best fitting model. The two fits coirespond to NGC 4550 (top) and NGC 4660 (bottom) respectively. 



in opposite directions. In all the fits we adopted a modest regular- 
ization (A = 10; see Ivan der Marel et alj jl998l) for a definition). 
In addition the SAURON data allow us to find that the two disks 
have very different scale-heights. Integrating the solution over the 
whole model we find that in NGC 4450 the counterrotating compo- 
nent constitutes about 50% of the total mass, while in the case of 
NGC 4473 the counterrotating component is about 30% of the total 
mass. A more detailed analysis of the incidence of kinematically 
distinct components in early-type galaxies goes beyond the scope 
of the present paper. 

For comparison with the two peculiar cases, in Fig.|4]we also 
show the best fit model for the protot ypical fast rotating disky ellip- 
tical NGC 4660 dBender et al.l(l994l) and the corresponding model 
solution in Fig. |5] We also show the model solution for the proto- 
typical slow rotating giant ellip tical NGC 4486 (M87; the bes t fit 
model was shown in Fig. 3 of 'Cappellari & McDermidll2005l) . It 
appears that the fast rotators tend to be characterized by flat stel- 
lar components (orbital starting conditions close to the equatorial 
plane) with angular momentum significantly smaller than the circu- 
lar orbit (starting conditions far from the bottom left/right corners of 
the plot). The relatively low angular momentum of these flat com- 
ponents produces a significant radial dispersion, which is the main 
reason of the observed anisotropy. The dominant radial dispersion 
in these flat disk-like stellar components is reminiscent of the one 
observed in disk galaxies, and generally attributed to disk heating 
processes (e.g. Shapi ro, Gerssen. & van der Mare l 2003). In some 
cases the disk-like components may happen to counterrotate, or two 
disks may be present, in which case tangential anisotropy dominates 
the observed anisotropy. 

Interestingly the disk-like components, which tend to charac- 
terise the fast rotators, not only seem to be distinct in integral space, 
but also differ in terms of their stellar population. As shown in 
Fig. 1 1 of iKuntschner et al., (.2006 , Paper VI), all the flat fast ro- 
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Figure 5. Integral space at a given energy for the solution of the 
Schwarzschild models for the galaxies NGC 4473, NGC 4550, NGC 4660 
and NGC 4486. Each panel plots the meridional plane (i?, z) with the loca- 
tion (white dots) where orbits are started with vj^ = Vz = al the given 
energy. For nearly edge-on galaxies, the white dots also corresponds to the 
position of the orbital cusps, where every orbit gives its strongest contribu- 
tion to the observables on the sky plane. The energy was chosen as that of a 
circular orbit with radius Rc = 32", which is about the size of the observed 
SAURON field (red rectangle). The coloured contours show the fraction of 
mass assigned to different orbits at that energy, where bright colours cor- 
responds to high mass fractions. Orbits at negative and positive R starting 
conditions correspond to prograde and retrograde respectively. Orbits with 
high angular momentum are found in the bottom right and bottom left cor- 
ner respectively on the plots. Orbits near the symmetry axis (green line) 
have low angular momentum. Orbits near the equatorial plane (z = 0) are 
intrinsically flat. See Fig. 6 of Paper IV for a detailed explanation of this 
diagram. 



tators are characterized by an Mgb line-strength distribution which 
is flatter than the isophotes. The fact that a flat Mgb is seen in all 
the flat galaxies suggests that perhaps all the fast rotators contain 
this metallicity-enhanced disk-like component, which is only vis- 
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Figure 6. Anisotropy parameter f3 in the meridional plane versus the 
luminosity-weighted second moment (Jc of the stellar velocity (from Pa- 
per rV). The red and blue labels refer to the slow and fast rotators respec- 
tively. 

ible in the flattest objects because th ey are closer to ed ge on. This 
result is reminiscent of the finding bv lLauer et al.l l l2005h that all flat 
cuspy galaxies show a disk. This metallicity enhancement indicates 
that additional star formation activity happened in the disk and, not 
surprisingly, implies that gas dissipation was involved in the forma- 
tion of these flat structures. The radial anisotropy however shows 
that heating was significant after the disk formation, as otherwise 
the stars in the flat components would still move on orbits that are 
closer to circular. 

4.6 Relation of anisotropy with other global observables 

In Fig. |6] we show the correlation between the anisotropy and the 
galaxy velocity dispersion CTc within the effective radius, taken from 
Paper IV. There is a trend for the most massive galaxies to have a 
smaller anisotropy within one Rc, with the exception of the special 
case NGC 4458 (see also Paper IV, § 5.1). In this diagram we use 
the parameter P to characterise the anisotropy. This parameter is 
measured in the galaxy meridional plane so it describes the orbital 
distribution in a way that is not affected by the direction of rota- 
tion of the stars in the galaxy. The (3 parameter measures the same 
anisotropy e.g. in two galaxies that formed with the same physical 
process, but in which one galaxy experienced a merger in a prograde 
direction and the other in a retrograde direction. 

A general trend is found between anisotropy and the intrin- 
sic galaxy ellipticity (Fig.|7}. The best-fitting relation to this rather 
scattered distribution of points, using a robust bisector algorithm, 
has the form: 

/3 = (0.6±0.1)ei„tr. (19) 

Given that the errors on both /3 and eintr are model dependent and 
difficult to estimate, we adopted constant errors on both variables. 
The quoted error on the slope was determined by enforcing the con- 
dition ~ ^, where u is the number of degrees of freedom of the 
fit. 

A better understanding of the meaning of the observed trend 
comes by considering the region in the (/3, eintr) plane that is al- 
lowed by the tensor virial equations. Equation l ll4t defines in fact a 
relation between the three variables eintr, V/a and 5. This equation 
is usually visualised in the observational plane {V/a, eintr), for dif- 
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Figure 7. Anisotropy parameter /3 in the meridional plane versus the average 
intrinsic ellipticity Eintr of the galaxies. The observed anisotropy is related 
to the intrinsic ellipticity of the galaxies. The thick black solid line is the 
best fitting relation. The dashed lines show the predictions of the tensor 
virial equations, in the form of equation )14t . for different values of the 
V/cT = 0, 0.25, . . . , 1 (see text for details). The red and blue labels refer 
to the slow and fast rotators respectively. The magenta line show the same 
relation /3 = 0.7eintr as the one in Fig.|9] 

ferent values of 5 as in Fig. [3] The same equation can be plotted in 
the theoretical plane [5, Eintr), for different values of V/a. As we 
saw that in real galaxies we have in general /3 ~ 5 (Fig.O, we 
will use the relation to understand the location of galaxies in the 
(/3, Eintr) plane (Fig.|7]l. As expected, within the numerical approx- 
imations and the limited spatial coverage, the galaxies generally lie 
within the region allowed by the tensor virial equations. The ob- 
served trend implies a specific form of the increase of the V/ a with 
increasing Eintr. 

A possible caveat is that a trend of the average anisotropy with 
increasing Eintr should be expected even if the galaxies were uni- 
formly distributed in f3 at every Eintr, due to the limits on P imposed 
by the virial equilibrium. However the 24 galaxies in the mod- 
elling subsample do not appear uniformly distributed in (3. More- 
over we show in the Appendix Icl that a uniform distribution in /3 
is not able to reproduce the observed distribution of all the galaxies 
of the SAURON sample on the {V/a, e) diagram: a general trend, 
and an upper limit of anisotropy as a function of ellipticity is re- 
quired to explain the data. As this trend is not implied by the ten- 
sor virial equations alone (nothing prevents all galaxies to scatter 
around /3 ~ 0), the observed distribution must hold a fossil record 
of the process that led to the galaxy formation. 

One explanation of the correlation may be that flatter galaxies 
were more strongly affected by dissipation before star formation 
was completed. Dissipation naturally makes a galaxy flat by reduc- 
ing the kinetic energy in the direction parallel to the global angular 
momentum. Dissipation is also required to produce the small ob- 
served kinematical misalignments for the fast rotators of the whole 
sample (Fig.[8}. 



5 {V/a, e) diagram OF THE SAURON SAMPLE 

In Section|4]we used dynamical modelling to recover the anisotropy 
of a sample of 24 galaxies extracted from the SAURON representa- 
tive sample. The galaxies in our sample were separated into two 
groups of fast and slow rotators (Paper IX). The slow rotators all 
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Figure 8. Kinemadcal misalignment ^! between tiie global photometric PA 
and the mean kinematic PA on the SAURON field. The red and blue labels 
refer to the slow and fast rotators respectively. We find that (i) all galaxies 
flatter than e ~ 0.2 have * ~ 0, and (ii) all the fast rotators have 'I' ~ 0, 
for any e (exception is the recent merger galaxy NGC 474). Assuming the 
SAURON survey constitutes a representative sample, this implies that triaxial 
galaxies have (c/a) > 0.8. Strongly prolate galaxies are also excluded from 
this sample. 



appear to have small anisotropy, while the anisotropy of the fast ro- 
tators seem to be related to the intrinsic ellipticity of each galaxy. 
One limitation of the analysis is due to the fact that the sample for 
which dynamical models were constructed is not necessarily repre- 
sentative of the whole galaxy population. To address this problem 
we show here that the findings of Section|4]can be extended to the 
whole SAURON sample of 48 E/SO galaxies. The results still holds 
even when we include in the analysis 18 additional E/SO galaxies 
that are not part of the survey, but for which SAURON data also 
exist. 



5.1 Intrinsic shape distribution 

A first point to address, to be able to simulate the observed distribu- 
tion of galaxies in the {V/a, e) diagram, is to understand what is the 
shape distribution of the galaxies in the SAURON sample. This is in 
principle a very difficult and intrinsically degenerate task, as even 
the full knowledge of the observed distribution of galaxy shapes 
and of their projected kinematical misalignments are not sufficient 
to recover the corres ponding intrinsic quantities, unless str ong as- 
sumptions are made jFranx. lUingworth. & de Zeeuwlll99lh . How- 
ever the availability of integral field stellar kinematics allows us to 
realise that the situation is simpler for real galaxies and that strong 
conclusions can be reached without any statistical analysis. 

In Fig. [8] we plot the kinematical misalignment ^ for the 
whole sample of 48 E/SO galaxies. ^ is defined as the angle be- 
tween the projected minor axis of the surface brightness and the 
direction of the projected kinematical rotation axis, which corre- 
sponds to the projected direction of the intrinsic angular momentum 
jFranx. lUingworth. & de Zeeuwlll99l]) : 



sin ^ = I sin(PAkin - PA. 



phot ) 



0° < * < 90° 



(20) 



The PAkin is defined as the PA at which the observed SAURON 
velocity field is best approximated by a bi-antisym metric version of 
it and was measured as described in Appendix C of lKrainovic et alj 
( 12006) . 



The position angle (PA) of the photometric major axis 
is obtained using the moments of the surface-brightness 
dCarter & Metcalfe! [19801) from the large scale MDM photometry 
(Falcon-Barroso et al. in preparation). It generally agrees well with 
the global PA quoted b y catalogues like R C3 and the 2MASS Ex- 
tended Source Catalog jjarrett et alj2000l) when available. It is im- 
portant to consider a global large-scale PA because: (i) early type 
galaxies tend to be r ounder in their central regions than in the 
outer parts (Rvde n. Forbes, & Terleviciill200lh . Thus the photomet- 
ric axes are generally more accurately determined at large radii, (ii) 
The central regions can be easily affected by bar perturbations. In 
the case of a disk galaxy with a bar we are not interested in the PA 
of the bar, but we want to measure the PA of the main di sk, which 
is better measured at large radii (e.g. lErwin et"aLll2003h . In most 
cases however, a constant-PA surface brightness distribution de- 
scribes well the whole galaxy and no distinction needs to be made. 

It turns out that nearly all the fast rotators in Fig. [8] have a 
misalignment ^' ~ 0° within the measurement errors. In the few 
cases where the misalignment reaches '3/ < 5° (e. g. NGC 3377, 
NGC 3384, NGC 4382, NGC 7332) the nonzero value can be as- 
cribed to the presence of a bar. The only exception is the galaxy 
NGC 474 for which ^ ~ 40° . This galaxy appears distorted by a 
recent merger, as evidenced by pro minent stellar shells at large radii 
jTumbull. Bridges. & Citfteill 19991) . Interestingly, also for the obvi- 
ous nearly face-on bars NGC 4477 and NGC 4262, characterized 
by the typical oval shape of the isophotes (see Paper III), the kine- 
matical PA appears not affected by the bar and is still well aligned 
with the outer disk, which defines the photometric PA. 

In a stationary triaxial system kinematical misalignments are 
common due to two reasons: (i) the intrinsic angular momentum 
can lie anywhere in the plane containing the long and short axis, so 
the kinematical axis are not necessarily aligned with the principal 
axes of the density; and (ii) projection effects can produce observed 
misalignments even when the angular momentum is aligned with a 
principal axis. The fact that ~ for all the fast rotators can then 
only be interpreted as strong evidence that these systems are nearly 
oblate systems. In some cases the fast rotators contain bars, which 
provide further evidence for their disk-like nature. 

The situation is completely different for the slow rotators. All 
these objects appear quite round in projection, the flattest one hav- 
ing e < 0.3 (excluding the special case NGC 4550). Some of these 
galaxies show significant kinematical misalignments, although in 
general the rotation is confined only to the central region, so it 
is not representative of the whole galaxy. But some galaxies like 
NGC 5198 can only be genuine triaxial systems (see Paper III). 
However the maximum observed ellipticity sets a limit of c/a > 0.7 
on the ratio between the shortest and longest axis of the density. A 
complete larger sample of slow rotating galaxies or detailed triaxial 
dynamical modelling would be needed to investigate the shape dis- 
tribution of these objects in more detail. The observed difference in 
shape between the fast rotating galaxies and the generally brighter 
slow rotating ones is consistent with the bimodality in the shape 
distribution between the fainter and brighter galaxies observed by 
Tremblav & Merritt ( 1996). The availability of the stellar kinemat- 
ics however shows that the shape distribution, not surprisingly, is 
more closely related to the kinematics than to the luminosity alone. 



5.2 Projection effects on the iV/a, e) diagram 

We established in Paper IX and further verified in Section [STI that 
the whole SAURON sample of E/SO galaxies can be broadly sub- 
divided into a class of fast rotating and nearly oblate systems and 
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another class of weakly triaxial slowly rotating galaxies. In Sec- 
tion |4j6] we saw that the anisotropy /3, measured in the meridional 
plane of the 24 galaxies in the modelling subsample, is roughly de- 
scribed by the relation (i ~ 0.6eintr. Moreover in Section 1421 we 
found that with a few notable exceptions [3^5, where 5 is the 
anisotropy as can be inferred from the {V/o, e) diagram. This im- 
plies that in most cases, one should also have 5 ~ 0.6 eintr. If these 
relations that we found for a small sample are valid in general, we 
should be able to model in a statistical way the distribution of the 
whole sample of fast rotators in the {V/a, e) diagram, as a random 
sample of oblate galaxies, with anisotropy defined by their intrinsic 
ellipticity eintr. The slow rotators should be inconsistent with this 
distribution, as they are not expected to be well described by oblate 
systems. A rigorous test of this idea is not possible due to the com- 
plex nature of our sample selection effects (Section [6.2t . however 
an attempt is still performed in Appendix |C] via Monte Carlo simu- 
lation. There we show that the observed distribution of galaxies on 
the iy/cT, e) diagram is indeed consistent with a nearly linear trend 
of anisotropy with ellipticity. 

Here we just try to obtain a qualitative understanding of the ef- 
fect of the inclination on the {V/a, e) diagram. For this we plot in 
Fig.|9]with a magenta lin^B a linear relatioi2|(5 = 0.7 eintr for edge- 
on galaxies. This relation approximately traces the lower envelope 
described by the location of the observed fast-rotating galaxies on 
the {V/o, e) diagram. It is steeper than the formally best-fitting re- 
lation l |19t . but still within the large errors. It is useful to have a 
qualitative understanding of the projection effects, as the paths fol- 
lowed by galaxies of different intrinsic (5, eintr) are nearly parallel, 
when the inclination is varied (e.g. Fig.O. For a given intrinsic el- 
lipticity Eintr and anisotropy 5, the corresponding V/a value, for 
an edge-on view is computed with the inverse of equation l ll4t . The 
dotted lines in Fig. |9] show the location of the galaxies, originally 
edge-on on the magenta line, when the inclination i is decreased. 
The projected V/a and ellipticity values at every inclination are 
computed with the inverse of equations l ll2l|13t . 

It appears that indeed all the fast-rotators lie on the {V/a, e) 
diagram to the left of the magenta line of the edge-on galaxies, as 
expected if they were a family of oblate models following a general 
trend of increasing 5 for larger e. The upper envelope of the region 
spanned by the lines of different inclinations is also generally con- 
sistent with the location of the observed galaxies, with reasonable 
values of the maximum intrinsic ellipticity and anisotropy. In gen- 
eral, the fast-rotators for which the inclination was determined in- 
dependently lie near a similar inclination on the diagram. The slow 
rotators on the other hand do not lie in the region allowed by oblate 
models following the given <5 — e relation. They may not follow 
any such relation or, more likely, they may not be oblate. This is 



^ A compact rational approximation of the magenta line which appears in 
the {V/(T, e) diagram of Fig.|9]is given by: 

which has a 1% accuracy in the useful interval e = [0.01, 0.7] . 

^ Any relation between anisotropy and ellipticity for oblate galaxies has to 

satisfy the necessary condition {V/a)'^ > in equation 1141 . or by Taylor 

expansion 

(5 < 1 - l/n{e) Ri 0.8£ + 0.15e^ + OMe^ + ... (22) 

This implies that the steepest linear relation is 5 = 0.8 e. independently of 
a. The magenta line is less steep than this limit. 




e 



Figure 9. (V/cr, e) diagram for the 48 galaxies in the SAURON representa- 
tive sample. The red and blue labels refer to the slow and fast rotators respec- 
tively (Table[T). The magenta line corresponds to the relation <5 = 0.7 Eintr 
for edge-on galaxies. The dotted lines show the location of galaxies, origi- 
nally on the magenta line, when the inclination is varied. Different lines are 
separated by steps of 10° in the inclination. The dashed lines are equally 
spaced in the intrinsic ellipticity. 



consistent with the result of Section [STI that the fast and slow ro- 
tating galaxies cannot be described by a single shape distribution. 
In Appendix Ic] we also discover that a sharp truncation below the 
magenta seems required to explain the observations. This zone of 
avoidance in the {V/a, e) diagram may be related to the stability 
of axisymmetric bodies with high ellipticity and large anisotropy, 
and needs to be confirmed with a larger and unbiased sample of 
galaxies. 

Remarkably, only two galaxies fall significantly below the 
magenta line for e > 0.3. These are the fast-rotator NGC 4473 
(marginally) and the special slow-rotator NGC 4550. This was ex- 
pected from the fact that these two galaxies are dominated by tan- 
gential anisotropy (Fig.|2j, due to the presence of two counterrotat- 
ing disks (Fig.|5]l. Although they follow the /3 — e relation, they do 
not satisfy the 5 — e relation as they have (3 < 5. 

Other significant exceptions are NGC 3156 and NGC 2685, 
which lie well above the general envelope described by the other 
galaxies, apparently next to the oblate isotropic line. Do these 
galaxies actually rotate as fast as an isotropic rotator? They are al- 
most completely disk-dominated, as evidenced by the nearly con- 
stant velocity dispersion maps, without clear signature of a hot cen- 
tral stellar component (Paper III). If this is the case they must be 
intrinsically quite flat and rotate significantly less than an isotropic 
rotator. For NGC 3156 the inclination i ~ 68° determined via 
modelling in Paper IV agrees with the inclination inferred from 
the shape of the outer isophotes, under the assumption of a thin 
disk geometry. At that inclination the models show that the galaxy 
indeed follows the general trend of anisotropy observed for the 
other galaxies (Fig. [TJ. A similar discussion is likely to apply to 
NGC 2685. These two galaxies appear to be special cases of flat 
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and disk-dominated galaxies which are poorly represented in our 
sample of early-type galaxies. 

The obs erved anisotropy of the fast-rotators is consistent with 
the finding bv B urkert & NaabI l l2005l) who recently pointed out that 
the observed distribution of disky ellipticals on the {y/o,e) dia- 
gram does not imply isotropy, but can be explained by projection 
effects on a family of anisotropic and nearly oblate merger rem- 
nants, as they obtained as a result of A'^-body simulations. 



5.3 E and SO galaxies on the {V/a, e) diagram 

In Paper IX we showed that early-type galaxies can be subdivided 
into two homogeneous classes of slow and fast-rotators respec- 
tively. The first class is composed almost exclusively of E galax- 
ies, while the second class is composed both of E and SO galaxies. 
Now SO galaxies are defined by the presence of a disk component 
and it is natural to ask whether all the fast rotators contain a disk, 
whic h is not always visib le photometrically due to projection ef- 
fects ( lRix&Whitj|l990h . In fact the disk becomes more difficult 
to detect when the galaxy is away from edge-on. To try to answer 
this question we can analyse the galaxies that, according to our sim- 
ple explanation of Section[5]2]are supposed to be close to edge-on. 
If they are indeed edge-on and if they all contain a disk compo- 
nent, this should be visible in the photometry. A detailed analysis 
of the photometric data will appear in a future paper in this series 
(Falcon-Barroso et al. in preparation) but we present here a qualita- 
tive discussion in relation to the {V/o, e) diagram. 

We find a range of properties in galaxies near the edge-on 
magenta line in Fig. |9] At the lowest V/o ~ 0.1 there is the 
SO galaxy NGC 4382, which shows evidence of interaction, stellar 
shells and distorted photometry at large radii from our MDM pho- 
tometry. At higher V/a ~ 0.2 we have NGC 821, NGC 4473 and 
NGC 4621, which are classified as E (Table [TJ, but show strongly 
disky isop hotes and would be classifie d as disky ellipticals E(d) ac- 
cording to lKormendv & Bendeij ( fl99db . An exception is NGC 2768, 
which does not have evidence for a disk. At higher V/a > 0.3 we 
find only SO galaxies. NGC 4270 has a peanut-shaped bulge, in- 
dicative of a nearly edge-on bar, NGC 3377 has strongly disky 
isophotes, while NGC 4570 and NGC 5308 have a clear thin disk 
component. 

Although the current galaxy sample is too small to reach firm 
conclusions, we speculate that the observed trend in galaxy prop- 
erties along the edge-on magenta line can be generalized to under- 
stand the whole distribution of the fast-rotators on the V/a dia- 
gram. The magenta line seems to describe a sequence of decreasing 
bulge over total-light ratio B /T. The region of low V/a is popu- 
lated by galaxies with B /T ~ 1. At increasing V/o the edge-on 
magenta line is dominated by SO galaxies, with B /T ~ 0.5 — 0.6 
fsimien & de Vaucoule urs 1986). When inclination effects are in- 
cluded on the magenta line one should expect the upper part of the 
envelope of the fast-rotators to be still dominated by SO galaxies, 
where the disk can still be detected even for a non edge-on view, 
while the lower part will be populated by non-disky E, where no 
evidence of the faint disk can be detected unless the galaxies are 
edge-on. This is what one can qualitatively infer from Fig. [TO] where 
we show the same sample as in Fig.|9] but with the morphological 
classification included. The decrease of B/T along the magenta 
line may also be the reason for the finding of Paper IX (their Fig. 7) 
that the fast-rotators with the largest (and thus generally the 
largest V/o) tend to be the least massive. 

The fact that the fast-rotators contain disks constitute a prob- 
lem for the characterization of their properties using the {V/a, e) 



b 

> 




Figure 10. Galaxy morphology on the {V/cr, e) diagram. Same as Fig.|9]but 
here filled triangles represent SO galaxies and open circles Es. The envelope 
of the fast-rotators is populated mainly by SO in the upper part and by E in 
the lower part. 



diagram. The diagram in fact is rigorously valid only for galaxies 
with a density stratified on similar oblate ellipsoids. This condition 
is clearly not satisfied for galaxies containing a rounder bulge and 
a flatter disk component. It would be interesting to understand in 
detail the orbital distribution of the different photometric compo- 
nents in early-type galaxies, as this constitutes a powerful record of 
the formation history. This would also be in principle feasible with 
the integral-field data and Schwarzschild models we have on a few 
carefully chosen galaxies, as the models do not make any assump- 
tion about the galaxy being made by a single stellar component (see 
e.g. Fig. Is}. This goes beyond the scope of the present paper. Here 
we notice that, from the dynamical models, there seems to be a ten- 
dency for the anisotropy to be stronger near the equatorial plane of 
the few disk-dominated galaxies (Section r4.lt . If this is true for the 
fast-rotators in general, one may speculate that the relation between 
anisotropy and ellipticity is related to the increase of both quantities 
with increasing B /T. However the anisotropy is in general not con- 
fined to the disks alone, indicating the situation is more complex. 

In this paper we do not attempt a disk-bulge decomposition 
as we focus on the statistical analysis of a large sample. We only 
consider the global properties of the stellar orbital distribution in the 
(possibly heterogeneous) region of the E and SO galaxies enclosed 
inside ~ li?e. Our goal is to study the ordered nature of early- type 
galaxies as seen on the (V/a, e) diagram, in much the same way 
as a similar set of simple measurements of integrated properties of 
a heterogeneous set of E and SO galaxies allows the Fundame ntal 
Plane to be studied (e.g. |j0rgensen. Franx. & Kiaergaardl 19961) . 



5.4 SAURON {V/o, e) diagram including 18 'specials' 

In addition to the 48 E/SO galaxies of the SAURON representative 
sample presented in Section [J!2l we observed 18 additional E/SO 
galaxies with SAURON, using the same setup as for the main survey. 
For these galaxies we extracted the stellar kinematics and measured 
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Table 3. Measured parameters for the 18 E/SO which are not part of the SAURON sample. 
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Notes: The meaning of the various columns is the same as in Table[T] Here the effective radii Rc are mostly taken from RC3. Exceptions are NGC 221 for 
which Rc was taken from Paper IV, and NGC 2679, NGC 6547 and NGC 6548, for which Rc was determined from our MDM photometry, 
t These two galaxies show nearly face-on bars. The ellipticity represents the one of the outer disk, which also defines the photometric PA. 



the {V/a, e) and Xr using the same procedures and definitions as 
for the main sample (Section |4.4t and give the numerical values 
in Table [3] In Fig. [TT|we show an extension of Fig. |9l where, in 
addition to the 48 galaxies of the SAURON survey, we also include 
the 18 specials. The new galaxies still populate the same regions as 
the galaxies from the main survey, confirming the result from that 
smaller sample. In particular all the fast-rotators still lie within the 
envelope defined by the S — eint relation (the magenta line) and its 
variation with inclination, while the slow-rotators lie in the bottom- 
left part of the ( V/cr, e) diagram. The five new slow-rotators are 
on average flatter than the ones in the survey sample, but still fall 
within the same ranges of low ellipticity (e < 0.3). 

A notable exception is NGC 720 which is at the edge of be- 
ing classified as a fast-rotator (see Table[3}, but is quite flat and lies 
well below the magenta line. This galaxy shows little rotation, but 
displays a clear and symmetric velocity pattern, with no kinemat- 
ical misalignment and no sign of photometric twist ( Peletier et al] 
[l990l) . All this indeed supports the classification of this galaxy 
as a fast-rotator and would suggest that the slow rotation may 
be due to counterrotating stellar components, as in the case of 
NGC 4473 and NGC 4550 (SectiongjSj. However, contrary to this 
hypothesis is the lack in NGC 720 of an elongation of the a field 
along the galaxy maj or axis. Two-integral dynamica l models using 
the Jeans equations l lBinnev, Davies. & Illingworth|[l99Q) already 
showed that part of the slow rotation in this galaxy must be in the 
form of a flattened velocity ellipsoid in the meridional plane [/3 > 
i.e. DF not of the form j(E, L^)]. This galaxy shows evidence from 
X-ra y observations of being embedded in a triaxial dar k-matter 
halo jRomanowskv & Kochane3ll998l : [Suote et al.ll2002h and the 
stellar component may be triaxial as well. Detailed three-integral 
Schwarzschild's dynamical modelling is needed to clarify the na- 
ture of the orbital distribution in this interesting galaxy. 

A key additional test provided by this extended sample con- 
cerns the shape distribution that we discussed in Section 15.11 In 
Fig. [12] we present a plot of the kinematical misalignment for the 




Figure 11. Same as Fig. [9] but for the 66 E/SO galaxies observed with 
SAURON (48 galaxies from the main survey plus 18 specials). The red and 
blue labels refer to the slow and fast rotators respectively. The filled sym- 
bols refer to the galaxies in the SAURON survey, while the NGC numbers 
indicate the 18 'special' objects. 



new sample of 18 galaxies, in addition to the SAURON sample. The 
sample of specials contain four additional slow rotators NGC 4168, 
NGC 4261, NGC 4365 and NGC 4406, which display a very clear 
rotation al ong the apparent minor axis and are ce rtainly not oblate 
bodies (see' Pavies et alj|200ll . lsta tler et al."2004' for a detailed dis- 
cussion about NGC 4365, but see lvan den Bosch et al..2007.) . These 
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Figure 12. Same as Fig.[8]but for the whole sample of 66 galaxies with exist- 
ing SAURON stellar kinematics. Blue and red symbols refer to the fast/slow- 
rotators respectively. The filled symbols refer to the galaxies in the SAURON 
survey, while the NGC numbers indicate the 18 'special' objects. The mor- 
phology of NGC 2679 is dominated by a strong face-on bar, which explain 
the observed large misalignment. All the other new fast-rotators are consis- 
tent with a family of nearly oblate bodies and show no significant kinemati- 
cal misalignment. 



new objects confirm that some of the slow-rotators must be tri- 
axial. A previous indication that large kinematics twists tend to 
appear only in the boxy s lowly rotating galaxies was found by 
iKormendv & Bended lll996h . The additional fast-rotators still show 
no significant kinematical misalignment (>]/ < 5°), the only excep- 
tion being NGC 2679 which has a morphology dominated by the 
obvious oval shape of a nearly face-on bar. Even the face-on barred 
galaxy NGC 6548 has a kinematical PA which well agrees with the 
PA of the large scale stellar disk. In total, out of 49 fast-rotators pre- 
sented in this paper, the only two that show significant kinematical 
misalignment, and are thus not consistent with axisymmetry, pos- 
sess strong bars or show signs of recent interaction. This provides 
strong evidence that all the fast-rotators constitute a homogeneous 
family of flattened and nearly oblate systems (sometimes with em- 
bedded bars). 



5.5 Caveats 

A number of caveats and possible improvements apply to the results 
we discussed in the previous sections: 

Sample size and selection: The sample of E/SO galaxies we con- 
sider in this paper is still relatively small and is not fully representa- 
tive of the whole galaxy population. The modest size of the sample 
constitutes by far the largest limitation of the present study. 

Field coverage: Although the available kinematics extends to 
larger distances and has higher S/N than what was used in pre- 
vious studies of the {V/o,e) diagram, it is still spatially limited 
and does not sample the whole galaxy, as would be required for 
a rigorous interpretation of the {V/o, e) diagram. The comparison 
between the modelling results, which treat the limited spatial cov- 
erage in detail, and the {V/(T,e) results (Section [4.4t . and the test 
with isotropic models (Appendix [BJ, all indicate however a general 
agreement, suggesting the results are reliable. Another effect is that 
slow-rotators tend to have larger Rc than the fast-rotators. For this 
reason the former are not spatially as well sampled as the latter. We 



verified that the conclusions would not change if we restricted our 
analysis to within Rc/2 for all galaxies. 

Bars: The SAURON representative sample includes a number of 
barred galaxies. Since their c omplex structure and dynamics (e.g. 
ISellwood & Wilkinson|[T993l) can not be represented in detail by 
oblate spheroids, and axisymmetric models likely provide only a 
very rough approximation, it is therefore unclear how bars affect 
correlations derived under the oblate assumption. Most strong bars 
must however have been rejected by the selection criteria for the 24 
galaxy in our modelling subsample, so they are unlikely to be the 
primary driver behind the observed trend between 13 and eintr. 

Triaxiality: The slow rotating galaxies are likely weakly triaxial 
(de Zeeuw & Franx 1991) and the assumption of axisymmetry in 
our models certainly introduces a bias in the recovered anisotropy. 
As it appears that the triaxiality of these systems is not strong, at 
least in their central regions, the inferred anisotropy is likely to be 
representative of what we would measure from a full triaxial model 
Jvan den Bosch et a l. 2007.) . A more detailed investigation of this 
issue will be presented in a future paper. 



6 COMPARISON WITH PREVIOUS WORK 

6.1 {V/ a, e) diagram from long-slit kinematics 

In this section we compare the results on the {V/a,e) diagram de- 
rived from our SAURON integral-field stellar kinematics to previous 
results obtained from long-slit kinematics. The (V/a, e) diagram 
fro m long- slit data of E g alaxies was fi rst presented and discussed 
by lUingworth ( 1977) and Binnev ( 1978). It was later extended to 



galaxy bulges by Kormend y & Illingworth ( 1982il. A widely know n 

Lay cm, 



Davies et ; 



form of the diagram is the one presented in 
which is based on a compilation of long-slit kinematical measure- 
ments for 50 elliptical gala xies and bu lges. The compilation was 
extended bv [Bender. Burstein. & F aber (1992) and many new ob- 
servations were added bv iBender et al. tl994.) . The largest homo- 
geneous set of long-slit V/a and ellipticity values is currently con- 
stituted by the 94 measurements for elliptical galaxies by R. Ben- 
der (private commu nication). We will refer to this sample as that of 
iBender et al. I ( Il994l) . although it was later extended with new un- 
published measurements. 

To perform a comparison between long-slit and integral-field 
V/ a measurements one has to take into account that the two tech- 
niques measure two different quantities and cannot be compared in 
a rigorous way. From integral-field data one measures a luminosity- 
weighted quantity (equation 1101 ), which is closely related to the 
corresponding quantity appearing in the tensor virial equations and 
can be used directly in a quantitative way. From long-slit kinemat- 
ics, one usually only measures the quantity Knax/o"o, where Vmax 
is the maximum observed velocity (usually on the galaxy major 
axis), while ctq is an average value of the velocity dispersion in- 
side Rc/2. This value cannot be linked to the expression appearing 
in the tensor virial equations without some rather arbitrary assump- 
tions or 'calibration' with models. 

From theoretical considerations, iBinnevI j2005h provides a 
very rough connection between the two V/a determinations. He 
shows that the long slit Knax/o"o measurements agree with the 
isotropic line, in the old formalism, if the isotropic line (his equa- 
tion [9]) is scaled by a factor 7r/4 so that for an oblate isotropic 
galaxy V/a = 0.99 at e = 0.5. In his new integral-field formalism 
the isotropic line is a factor ^/2 lower (his equation [26]) than that 
in the old one when a = 0, so that for an oblate isotropic galaxy 
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Figure 13. Comparison between the SAURON integral-field (V/(j)c mea- 
surements and the corresponding Vmax/co values by Bender et al. 1 1994}), 
corrected with equation (23) to make them comparable to integral-field ob- 
servations. 



V/a = 0.89 at £ = 0.5. For our adopted value a = 0.15 (Ap- 
pendix |B]i one has V/a — 0.79 at e = 0.5 so that the expected 
combined decrease from long-slit data to integral-field one would 
be a f actor (V / a)e ~ 0.8 (Vmax/cro). However, as also pointed 
out by iBinnevI j2005l) . the precise form of this equation depends 
on the detailed velocity distribution in a galaxy, as well as on the 
field coverage of the observations, and should not be expected to be 
quantitatively accurate. This was in fact his motivation to introduce 
a new formalism for integral-field data. 

To avoid the approximations and assumptions in the above the- 
oretical derivation, and to be able to compare directly our SAURON 
V/a values with th ose measured with long-slit observations by 
iBender et al. I ( Il994h . we selected the 28 galaxies in common be- 
tween the two samples and we determined the best scaling factor 
by least-squares fitting (Fig.|13l>. We found a best-fitting relation 



'Y.) R,0.57— 

\ a J a ao 



(23) 



which is ~ 30% lower than the above rough theoretical estimation. 
Although the differences in individual cases are not insignificant, 
there appears to be a general agreement between the SAURON and 
(scaled) long-slit V/a determinations. The main exception are the 
galaxies with V/a ^ 0.1, where the long-slit data tend to measure 
systematically smaller valu es. A comparison f or the ellipticity val- 
ues measured by us and by Bender et al.l ( ll994l) is shown in Fig.ll4l 
There is a very good agreement between the ellipticities measured 
with our different methods. 

In Fig.[T5]we plot the 94 long-slit V/a measurements by Ben- 
der, rescaled as in equation l l23t , and his corresponding ellipticity 
determinations, on top of the same lines as in Fig. |9] For the 28 
galaxies in the Bender sample for which we have SAURON obser- 
vations, we also show with an arrow the location of the correspond- 
ing SAURON V/a and e values. This plot shows that the use of 
integral-field data does not systematically change the distribution 
of galaxies on the (V/a, e) diagram. The differences between the 
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ity £ derived for th e present paper and the coiTesponding £394 values by 
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SAURON {V/a, e) diagram and the one in Fig.[T5]are discussed in 
the following section. 

The comparison of this section shows that the long-slit obser- 
vations already provided a reasonable approximation to the more 
rigorous V/a that can be obtained from integral-field kinemat- 
ics. However integral-field kinematics provides many advantages 
over long-slit measurements for the study of the {V/a, e) diagram: 
(i) The dynamical models can only be uniquely constrained by 
integral-field kinematics; (ii) V/a can be compared to the theoreti- 
cal value in a rigorous way, without the need for an ad-hoc scaling, 
allowing for a quantitative comparison with modelling results, in in- 
dividual cases, as we did in Section l4!4l (iii) the ability to accurately 
determine the kinematical misalignment allows us to recognise the 
different shape distribution of the fast/slow-rotators; (iv) the possi- 
bility to rigorously measure the specific angular momentum adds 
the crucial ability to distinguish between fast and slow-rotators, in- 
dependently from the {V/a, e) diagram itself (Paper IX). A proper 
kinematical classification of the galaxies is the key for the interpre- 
tation of the {V/a, e) diagram. 



6.2 Sample selection effects on the {V/a, e) diagram 

The comparison between the SAURON {V/a,e) diagram of Fig. [9] 
and the corresponding one bv lBenderet"ai]jl994h in Fig.|I5|shows 
clear differences in the distribution of the observed galaxies. The 
reason for these differences are easily understood as due to very 
different selection criteria in the two samples. 

The early-type galaxies of the SAURON survey were selected 
to have an absolute blue magnitude Mb ^ — 18, and in this lu- 
minosity range to include 24 objects in both the E and SO mor- 
phological classes, to be unif ormly sampled in Mb and in the ob- 
served ellipticity (as given bv lde Vaucouleurs et alJ[l99Tl . hereafter 
RC3). By construction, this sample provides a uniform coverage 
of the {V/a,e) diagram, but does not give a representative sam- 
ple of early-type galaxies in the given magnitude range. In fact the 
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Figure 15. (V/cr, e) diagram from long-slit kinematics. Same as Fig. |9] 
but for the 94 Vma x /p-q long-slit measurements and e determinations by 
iBender et alj jl994'). The Vmax/cro values have been corrected with equa- 
tion {23) to make them comparable to integral-field observations. For the 28 
galaxies for which SAURON observations also exist, the arrows indicate the 
location on the diagram of the corresponding value of Fig. [171 The red and 
blue lines refer to the slow and fast rotators respectively, classified according 
to the SAURON kinematics (Table[T). 



true observed shape distribution of E/SO galaxies displays a clear 
maximum at an ellipticity e ~ 0.25 and lacks very ro und and very 
flat galaxies fe.g. lLamhas, Maddox. & Lovedavll992h . The require- 
ment of a uniform sampling in e obviously produces an excess of 
very round and very flat galaxies with respect to the general galaxy 
population. An additional complication arises from the fact that 
the SAURON selection was done according to the RC3 ellipticities, 
which are measured at large radii. The ellipticity in galaxies is gen- 
erally not constant and its variation depends on galaxy luminosity, 
in the sense that luminous galaxies ten d to become rounder in their 
central parts than in their outer parts /Rvd en. Forbes. & TerlevichI 
liooi). The e used in the (V/a, e) diagram are those measured in 
the central regions, where the kinematics is available, and not those 
of RC3. As a result of this, the SAURON selection produces a distri- 
bution in e which is difficult to quantify. 

Although the SAURON sample selection, by construction, 
erases any information given by the e distribution, one may naively 
expect the distribution in V/a for any given e interval, to be rep- 
resentative of the general E/SO population, as no explicit selection 
was performed in V/o. In practice, this is not the case because of 
the morphological selection, which implicitly introduces a selec- 
tion on V/a. In particular Table[T]shows that out of the 24 galaxies 
classified as E from RC3, 14 (58%) are classified as fast-rotators 
and 10 (42%) as slow-rotators. For the SO galaxies 23 (96%) are 
fast-rotators and only I (4%) is a slow-rotator. Fig.|9]shows that the 
fast-rotators tend to lie above the magenta line, at any given e, while 
the slow-rotators are in the bottom-left part of the iV/cr, e) diagram. 
This implies that the selection of SO galaxies will guarantee objects 
above the magenta line in the {V/o, e) diagram, while E galaxies 
will provide an equal number of objects above the magenta line and 



in the bottom-left part of the diagram. The morphological selection 
thus defines the allowed ranges of V/ a of any galaxy. The require- 
ment of an equal number of E and SO galaxies in every interval of e 
specifies a certain distribution of galaxies in the (V/a, e) diagram, 
which is not representative of the early-type galaxies population. 
These selection biases are very difficult to estimate and to correct 
for from the present sample, as they depend on the unknown distri- 
bution of fast/slow-rotators as a function of intrinsic galaxy shape 
and morphological classifications. 

The morphological selection also affects the appearance of the 
(V/a,e) diagram a s obtained from the long-slit measurements by 
iBender et alT ( ll994i) . That sample is in fact constituted only by E 
galaxies. As originally pointed out by Davies et al. ( 1983), clearly 
stated by iKormendv & Bendeil jl996h and strongly confirmed in 
this paper and Paper IX, the E class constitutes a quite heteroge- 
neous group, made by less than half of genuine non-rotating ellip- 
ticals (the slow-rotators) and the rest by disk-like systems, resem- 
bling SO galaxies (the fast-rotators). The appearance of the (V/a, e) 
diagram of E galaxies thus depends on the fraction of fast-rotators 
that are (mis)classified as E in different e intervals. A direct com- 
parison with the SAURON (V/a, e) diagram for the E galaxies alone 
(open circles) can be seen in Fig.[TO] Although the SAURON sam- 
ple has only 25% of the number of E of .Bender et aL CI994.) , the 
distribution of the galaxies appears qualitative ly similar. 

The (V/a, e) values for the E galaxies of lBender et alj ( Il994h 
in Fig.[T5]span the same region as those of SAURON of Fig.l9lorlIll 
However there are three main qualitative differences between this 
diagram and the SAURON one, which can be understood as due to 
differences in the sample selection: (i) the E selection in the long- 
slit sample poses an upper limit on e < 0.5, as otherwise the galax- 
ies would be classified as SO. (ii) The lack of SO in the long-slit 
sample implies a smaller fraction of fast-rotating galaxies (above 
the magenta line), compared to the SAURON {V/a, e) diagram, (iii) 
The fact that there was no explicit selection on e in the long-slit 
sample gives a smaller fraction of very round galaxies, compared to 
the SAURON sample. Moreover older samples of V/a values tend 
to prefer flatter E, as these were thought to be most likely to be seen 
edge-on, and consequently better consistent with the assumption of 
the {V/a,e) diagram. 

It is important to obtain an unbiased view of the {V/a, e) dia- 
gram, from a larger statistical sample, to be able to compare the ob- 
served distribution of fast-slow rotators, with the predictions from 
different galaxy formation scenarios. For this, one needs a sam- 
ple which is only selected in absolute magnitude and is not bi- 
ased by the morphological selection, as the currently available sam- 
ples. Some selection in morphological type is still required, as the 
{V/a, e) diagram looses its meaning for spiral galaxies, which have 
a clear bulge and a thin disk component with very different elliptic- 
ity. However, contrary to the problematic classification into E or SO 
galaxies, the separation bet ween early-types galaxies and spirals i s 
robust and clearly bimodal tStrateva et alj200ll : lBaldrv et alj2004l) . 
Only with this larger and unbiased sample of E/SO galaxies will we 
be able to understand whether the fast-rotators all really fall within 
the envelope indicated by the current study, and whether the distri- 
bution of the slow/fast-rotators shows a bimodality in the (V/a, e) 
diagram, as the current data seem to suggest. 

6.3 Spurious increase in anisotropy of slow-rotators? 

In Section 16.21 we showed that our SAURON data on the 
{V/a,e) diagram are generally consistent with previous long-slit 
works. The current study however seems inconsistent with pre- 
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vious studies (e.g. iDavies et al ] |1983I: Ib innev &Tremainelll98l 
iKormendv & Bended 1 199^ in its finding that the massive and 
slowly rotating galaxies are not more anisotropic than the generally 
fainter and fast-rotating galaxies (e.g. Fig.|6j. 

Part of the difference between the current work and the pre- 
vious ones may lie in the sample selection, which here includes 
SO galaxies, as they cannot be distinguished from the other fast- 
rotators. The rest can be understood as due to (i) the use of the 
parameter (V/a)* as a measure of galaxy anisotropy and (ii) the 
neglect of inclination effects on the {V/a, e) d iagram. 

The anisotropy parameter was defined byl&rmen^ |l982B); 
iDavies et al.l ( ll983l) ; lBender & Nietol ( Il99(]h as 



(V/a) 



V/a 



(V/c 



(24) 



i.e. the ratio between the observed V/a value for a galaxy and the 
predicted value (V^/a)iso for an oblate galaxy with isotropic ve- 
locity dispersion tensor (Binnev 1978). As pointed out by Binne^ 
f200i . (V/a)* is not a reliable way of measuring anisotropy in a 
galaxy and should be replaced by the real anisotropy value S. The 
problems of (V/a)* as a measure of anisotropy S are the follow- 
ing: (i) galaxies with the same 5 but with different intrinsic ellip- 
ticity Eint, can appear to have dramatically different (V/a)* when 
seen edge-on, and even more different log{V/a)*. In particular the 
rounder galaxies will have a systematically lower (V/a)* at any 
given S (Fig. I16t. (ii) A galaxy with a given 5 and eint, will ap- 
pear to have a lower (V/a)* at lower inclinatio ns (i = 90° be ing 
edge-on) as shown e.g. in Fig.[9](see also Burkert & Naabll2005l) . 

The slow-rotators generally have observed ellipticities e < 0.3 
(Fig. |9l and 1151), which implies they are intrinsically not very flat. 
The fast rotators can be intrinsically flat. Given this difference in 
the shape of the fast and slow-rotators, both the ellipticity and in- 
clination effects mentioned in the previous paragraph act to sys- 
tematically decrease the observed (V/a)* of the slow-rotators with 
respect to the fast-rotators, giving the misleading impression of a 
strong increase of the anisotropy. We conclude that (V/a)* should 
not be used as a measure of anisotropy. 



7 DISCUSSION 

7.1 Early-type galaxies in the nearby universe 

Observations at high redshift, combined with detailed measure- 
ments of the cosmic microwave background, have provided a sce- 
nario for galaxy formation in which the Universe is dominated by 
dark matter of unknown nature. The general picture for the assem- 
bly of the dark matter seems to be rather well understood and is 
described by a hierarchical process driven by the sole influence of 
gravity (e.g. ISpringel et aLii2005.) . Much however still has to be 
learned about the way the luminous matter, namely the stars and 
galaxies that we actually observe, form from the accretion of gas 
within the dark matter potential. The complication comes from the 
fact that these processes are driven by the hydrodynamic and chem- 
istry of the gas, combined with complex feedback processes of the 
radiation on it. 

One way to constrain the formation of the luminous matter 
in galaxies is by trying to uncover the 'fossil record' of that pro- 
cess in nearby galaxies, which can be observed in great detail due 
to their vicinity. A key element to differentiate between different 
galaxy-formation model is the degree of anisotropy and the amount 
of rotation in early-type galaxies, in relation to other galaxy prop- 
erties like luminosity, surface brightness profiles or chemical com- 





> 
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Figure 16. Top Panel: (V/cr, e) diagram. The meaning of the solid lines is as 
in Fig.|9] The blue diamonds and red squares represent the location of edge- 
on galaxies, with different intrinsic ellipticity e^^i but the same anisotropy 
5 = 0.1. The red ellipse indicates the region where the slow-rotators are 
generally observed (Fig.|9land |15) . Bottom Panel: Location of the galaxies 
of the top panel in the (log(y/ a)* ,e) diagram. The diamonds and squares 
represent galaxies with the same small anisotropy, so one should expect 
\og(V /cr)* Ki 0. There is however a dramatic decrease oflogiV/cr)* for 
the objects which fall close to the region of the slow-rotators (red squares), 
which is not associated to any variation in the anisotropy. 



position (e.g. lBinnevlll978l : IDavies et al.lll983l) . Observations show 
that early-type galaxies can be broadly separated into the most mas- 
sive giant ellipticals on one side, which are red, metal-rich, slowly- 
rotating, and have shallow central photometric profiles, and the 
fainter objects on the other side, which tend to be bluer, metal-poor, 
are do minated by rotation, have cuspy profiles and may contain 
disks toavies et al. "l983'; 'Jaffe et al.* fl994l: iKormendv & Bendej 



1996; Faber et al. 1997; Lauer et al. 2005 

A first factor determining these observed differences in the 
structure of early-type galaxies appears to be the amount of dissipa- 
tion during the last galaxy merger. In particular more anisotropic 
and slowly rotating galaxies result from predominantly coUi- 
sionless major mergers, while faster rotating galaxies are pro- 
duced by mo re gas-rich merge r s, where dissipation plays an im- 
portant rol e jKormendvl 1 1 9891 ; iBender, Burstein, & Fabe^ 1 19921 ; 
iFaberet alj|l997l) . The prevalence of g as-poor mergers in the hier - 
archical formation of massive galaxies l lKhochfar & Burkerlll200^ 
would help explain the fact that slowly rotating galaxies are more 
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common in the high-mass galaxy range. A second factor driv- 
ing the galaxy structure seems to be the mass fraction of the 
merger components. Equal-mass galaxy mergers tend to produce 
more slowly ro tating galaxies than mergers of significantlv differ- 
ent m ass ratio jNaab. Burkert. & Hernquis^l 19991 ; iNaab & BurkertI 
I2OO3I) . A third factor appears to be the feedback on the gas, 
during a merger event, produced by the radiation emitte d by a 
central AGN or by supernovae winds (e.g. Granato et al.l |2004| : 
iSpringel. Pi Matteo. & Hernguis tll2005l) . 

A complete understanding of the formation of early-type 
galaxies, and the relative role of the different factors, requires the 
treatment of the above phenomena in a cosmological context, within 
the hierarchical merging process. This still has to be attempted in a 
rigorous way as current computer simulations cannot resolve the 
many orders of magnitude scales that are involved in these pro- 
cesses. This implies that instead of describing the gas via accurate 
hydrodynamical and radiative transfer equations, the models use ap- 
proximate and very simple relations, based on empirical relations, 
to describe the process. This paper and the companion Paper IX 
provide new observational constraints for the upcoming models. In 
the next section we speculate on the formation scenarios for the ob- 
served characteristics of early-type galaxies that we derived. 



7.2 Nature of the fast and slow rotators 

The fast rotators seem consistent with galaxies with a significant 
disk component, which experi enced minor mergers and accreted 
a significant amount of gas /Naab, Burkert, & Hernquistj 1 19991: 
Naab & Burkert 2003; Bournaud, Jog, & Combes 2005; Co x et alj 
20061 ; IJesseit et alj 120071 ; LNovak et al. 2006). The accreted gas. 



which settles in the galaxy equatorial plane and forms stars, would 
explain the metallicity enhancement in the plane of the disk (Pa- 
per VI) and the nearly perfect alignment of the angular momen- 
tum of the different, sometimes counterrotating, stellar structures 
(Section|43}. The minor mergers may explain the fact that rounder 
galaxies have a larger dispersion orthogonal to the disk plane 
(along the z-axis). The disk instabilities triggered by minor mergers 
can in fact deviate stellar orbits out of the disk plane and produce an 
increase in the Gz dispersion and a corresponding thickening of the 
disk. The trend may also be related to the two-components nature 
of the fast-rotators. Disks may be more anisotropic than bulges and 
are more prominent in flat galaxies (Section [53}. Numerical sim- 
ulations seem able to produce the observed level of anisotropy in 
these systems dBurkert & Naab 2005), but it is unclear for what con- 
ditions they will reproduce the observed trend between anisotropy 
and flattening. 

Although we do not find any physical distinction between E 
fast-rotators with disky isophotes and SO fast-rotators with an ob- 
vious disk component, and this suggests a common formation sce- 
nario for the whole class of objects, the current observations do not 
necessarily require a common origin. 

Originally the slow rotators were thought to form at hig h red- 
shift via equal mass mergers of disk galaxies jBamesiri988h . The 
N-body simulations were able to produce some slowly rotating sys- 
tems that appeared like ellipticals, however it became clear that this 
simple picture could not explain the observed appearance of slowly 
rotating ellipticals as a class (Hernquist 1992; Naab & Burkert 
l2003h . In addition, current numerical simulations of gas-free ma- 
jor mergers seem to generate large fractions of strongly triax- 
ial or prolate galaxies (0.2 < e < 0.6) than we do not observe 
(e.g. Bovlan-Kolchin, Ma, & Ouataert 2005;'Burkert & Naa bl2()0^ ; 
Icionzalez-Garcfa & van Albada.2005. ; .Cox et al.,,2006,) . A missing 



ingredient in these models is the cosmological context, which pro- 
duces sequences of mergers. Also lacking is a detailed treatment of 
the gaseous component, which is likely to play an important role in 
the early universe. Major gas-rich mergers, regulated by the feed- 
back of a powerful central AGN, have b een suggested to explai n the 
scaling relations of early-type gal axies teobertson et alj|2006l) and 
to rep roduce the dichotomy (e.g. 'pabe r et alj|2007l ; ICattaneo et aj] 
I2OO6I) in the colour distri butio n of g alaxies as a f unction of their 
luminosity ( Strateva et al.ll200ll ; lBaldrv et alj|2004f) . 

It is unclear why the slow rotators in our sample possess the 
low degree of anisotropy that we observe in their central regions. 
A process that could in principle explain both the relatively small 
anisotropy and the roundish shape of the slow rotators would be 
the global transition t o chaos, as can be produ ced by central super- 
massive black holes jMerritt & Fri dmanl4996). Current estimates 
of black holes mass fractions iHaring & Rix 2004) seem however 
too small to produce significant effects. The transition to chaos and 
the symmetrisation of the major triaxial merger remnant could be 
significantly enhanced by the gas inflow and its nuclear concen- 
tration during the initial phase of the mergers, be fore the gas is 
expelled from the galaxies by the AGN activity l lBamesl ll988). 
Another process that could explain the isotropy and shape of the 
slow-rotators, especially for the most massive galaxies found near 
the centre of cluster potential s, is a sequence of radial mergers 
jCoorav & Milosavlievidl2005h . 



7.3 Connection with higher redshift observations 

The detailed observations of the structure of the small sample of 
nearby galaxies that we discuss can be related to much less detailed 
observations at higher redshift obtained with much larger survey. 
From the analysis of a large sample of galaxies in the nearby uni- 
verse, observed by the Sloan Digital Sky Survey, it was found that 
galaxies are well distinct in terms of their colour, in the sense that 
they can b e clearly separated in a so called 'blue c loud' and a 'red 
sequence' jStrateva et alj200ll ; lBaIdrv et al.ll2004h . This discovery, 
and the subsequent realisation that the bimodality can be traced 
back in time at higher and higher redshift teell et al.l20"04h . allowed 
a dramatic improvement in the detailed testing of galaxy forma- 
tion scenarios. It is now believed that the gas-free and dead, red- 
sequence galaxies form by merging of gas-rich, star-forming, blue- 
cloud galaxies, followed by a rapid ejection of the gas , due to the 
feedb ack from a central supermassive black hole (e.g. Ipaber et ^ 
I2OO7I) . In a merger event between blue gas-rich galaxies the gas 
tends to settle onto a plane and form a disk, so that the end result 
of the merger, after the gas has been expelled from the system, will 
be a red, purely stellar system dominated by rotation. But mergers 
can still happen between red gas-poor galaxies, in which case sim- 
ulations show that the resulting red galaxy will show little or no 
rotation. 

We speculate that the two classes of fast and slow-rotating 
galaxies that we described may be the relics of the different forma- 
tion paths followed by early-type galaxies when loosing their gas 
content and moving from the 'blue cloud' to the 'red sequence', 
in the 'quenching' scenario of lFaber et all j2007h . In that picture 
the fast-rotators are naturally associated with the spiral galaxies 
whose disks were quenched by ram-pressure stripping or other ram- 
starvation processes, and subsequently experienced mainly minor 
mergers. The slow-rotators can be generated by a sequence of merg- 
ers involving little amount of gas along the red sequence, or due to 
major violent gas rich mergers in which the gas component was 
rapidly expelled by a starburst or a central AGN. 
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8 CONCLUSIONS 

We analysed the orbital distribution of elliptical (E) and lenticu- 
lar (SO) galaxies using observations of the stellar kinematics within 
one effective (half light) radius obtained with the SAURON integral- 
field spectrograph. We constructed the (V/a, e) diagram for the 48 
E/SO of the SAURON sample. For a subsample of 24 of the galaxies, 
consistent with axisymmetry. we use three-integral axisymmetric 
Schwarzschild models to recover the detailed orbital distribution. 

We used a new classification of early-type galaxies based on 
the stellar kinematics (Paper IX), to divide the sample into two 
classes of objects with and without a significant amount of angu- 
lar momentum per unit stellar mass. We refer to these two types of 
object as "slow-rotators" and "fast-rotators" respectively. The slow- 
rotators are more common among the brightest systems, and are 
generally classified as E from photometry alone. The fast-rotators 
are generally fainter and are classified either E or SO. 

From the combined analysis of the dynamical models and 
the {V/a, e) diagram we find that, in the central regions (within 
~ 1-Rc) of the galaxies which we sample with our kinematics: 

(i) The slow rotators can display misaligiraient between the pho- 
tometric and kinematic axes, indicating that as a class they are triax- 
ial. Those in our sample tend to be fairly round (e < 0.3), and span a 
range of moderate anisotropics (Binney 1978 anisotropy parameter 
5<0.3). 

(ii) The fast-rotators can appear flattened (e < 0.7), do not show 
significant kinematical misalignment (unless barred or interacting), 
which indicates they are all nearly axisymmetric, and span a large 
range of anisotropics ((5<0.5). The main physical difference be- 
tween fast and slow-rotators is not their degree of anisotropy, hut 
their intrinsic shape. 

(iii) The measured anisotropy appears to be due to a flattening 
of the velocity ellipsoid in the meridional plane {gr > Oz), which 
we quantify with the P anisotropy parameter. While the velocity 
dispersion tends to be similar along any direction orthogonal to the 
symmetry axis {an ~ (t^). 

(iv) We find a general trend between the anisotropy in the merid- 
ional plane j3 and the galaxy intrinsic ellipticity eintr. Flatter galax- 
ies tend to be more anisotropic. 

(v) The results from the models are consistent with the values 
of the anisotropy inferred from the ( V/cr, e) diagram of the whole 
SAURON sample of 48 E/SO, which we determined from our inte- 
gral field kinematics. We find that the distribution of the fast rotat- 
ing galaxies is well described as a family of nearly oblate galaxies 
following a trend between flattening and anisotropy. This result is 
confirmed when we include in the sample 1 8 additional E/SO galax- 
ies observed with SAURON. 



Although the sample of early-type galaxies considered in this 
work constitutes the first significant sample for which the {V/cy, e) 
diagram and the anisotropy could be investigated with integral-field 
stellar kinematical data, the main limitation of our work is the rel- 
atively small number of galaxies and the complicated sample se- 
lection effects. Even so, as we discussed in Section|7] we are able 
to reach some provisional conclusion on the formation of early-type 
galaxies. We reviewed the usefulness of the {y/o,e) diagram to test 
galaxy-formation models. A similar analysis, performed on a much 
larger sample of galaxies could provide much stronger constraints 
on the process by which galaxies form. 
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APPENDIX A: TESTING ANISOTROPY WITH 
TWO-INTEGRAL MODELS 

The three-integral axisymmetric models of Section 143] show that, 
in the central regions that we constrain with our integral-field kine- 
matics, (i) the flatter fast-rotators are anisotropic, in the sense that 
their velocity ellipsoid is flatter along the direction of the symme- 
try axis than along the radial direction < ctr; (ii) the rounder 
slow-rotators are nearly isotropic. In this section we do a qualita- 
tive test of the consistency and robustness of these statements with 
the results of two-integral Jeans models. This test uses a very dif- 
ferent and simpler modeling method and a different kinematics ex- 
traction. It allows us to verify the robustness of our results against 
implementation and observational details. 

The Jeans models we use assume axisymmetry, constant M/ L 
ratio, and a stellar distribution function DF of the form / = 
f{E,Lz), where E is the energy and Lz is the angular momen- 
tum along the symmetry axis z. From the same MGE parameter- 
ization of the surface brightness used in the Schwarzschild mod- 
els (Section [3j and an estimate of the central black-hole mass, 
we compute the unique predictions for the velocity second mo- 
ments /J2, at the best -fitting inclination, using the formalism of 
lEmsellem et iD (Ll994). See Paper IV for more details. The second 
velocity moments were extracted from the SAURON kinematics as 
jj,2 = VV^ + o^, where V and a are the mean velocity and the 
velocity dispersion used also in the (V/a, e) diagram (Section[2](. 

Al Anisotropy of the flat fast-rotators 

To test the anisotropy of the flat fast rotators we consider the three 
galaxies NGC 821, NGC 3377 and NGC 4660, which are con- 
strained by the observed photometry to be close to edge- on. This 
choic e reduces the non-uniquness in the mass deprojection jRvbickil 
Il987[) and correspondingly the possible degeneracies in the depro- 
jection of the higher velocity moments. These three galaxies have 
a a field well resolved by the SAURON instrumental resolution. For 
edge-on galaxies one does not have any direct measurement of the 
component of the velocity dispersion CTz parallel to the symmetry 
axis, however the flattening of the velocity ellipsoid is tightly con- 
strained if the galaxy is assumed to be axisymmetric and in equi- 
librium (same as in the {V/a, e) diagram). Nearly edge-on galaxies 
are thus ideal for our test. 

The comparison with the Jeans model predictions is shown in 
Fig. lAll The models provide a reasonable first-order description of 
the observations, however significant systematic deviations exists. 
This indicates that, assuming the mass traces the light, the modelled 
galaxies do not have a DF of the form f{E, Lz). In particular a 
common characteristics of models is that for all three galaxies they 
predict too much motion (high /12 values) along the galaxies pro- 
jected major axis with respect to the minor axis. Any two-integral 
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Figure Al. Jeans models of flattened galaxies. Left panels: bi-symmetrised 
and linearly interpolated maps of the observed second velocity moment 
112 = \/V'^ + cr^, measured from the SAURON stellar kinematics of the 
galaxies NGC 821. NGC 3377 and NGC 4660. The actual Voronoi bins 
in which the kinematics was measured ai'e shown in Paper III. All three 
galaxies are quite fattened and are constrained by the observed photometry 
to be close to edge-on. Right panels: two-integral Jeans modelling of the 
unique projected second moments associated to the MGE mass density dis- 
tribution, at the best fitting inclination from Paper IV. In all three cases the 
models show a larger difference between the /i2 measured along the major 
and minor axes than is observed. 

,f{E, Lz) model has a velocity-ellipsoid which appears circular in 
cross section with the meridional plane (aR = (Jz and {vzVr) = 0). 
The observed difference between models and observations indicates 
that the real galaxies have a larger velocity dispersion than the mod- 
els, along the R direction, orthogonally to the z symmetry axis. 
This implies that the velocity ellipsoid must be on average elon- 
gated in the 7? direction. In this way, for nearly edge-on galaxies, 
the observer sees a larger component of the dispersion along the 
minor axis, when the radial direction in the galaxy is parallel to the 
line-of-sight , than along th e major axis. A similar finding was in- 
terpreted by van der Marell ( [1991) as evidence of radial anisotropy 
and the same reasoning applies in our case. 

A2 Near isotropy of the round slow-rotators 

To test the near isotropy inferred using three-integral models for the 
rounder slow-rotators of the modelling subsample, we selected the 
three apparently quite round galaxies NGC 4374, NGC 4486 and 
NGC 4552. These three galaxies have none or very little evidence 
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Figure A2. Jeans models of nearly round galaxies. The open squares 
with error bars represent the observed second velocity moment /i2 = 
\/V'^ + cr'^, measured from the SAURON stellar kinematics of the galax- 
ies NGC 4374, NGC 4486 and NGC 4552. The values and errors have been 
computed as the biweight mean and standard deviation of the kinematics 
extracted from the Voronoi bins within circular annuli. The solid line shows 
the second velocity moment computed from two-integral Jeans models. The 
dashed vertical line shows the region where the SAURON field provides a 
full spatial coverage of the galaxy. This is the region visualized in Fig. [T] 
and from which the mean anisotropy was estimated. The two-integral mod- 
els provide a very good description of the data within the central region. 
Given the near spherical symmetry the models are almost fully isotropic in 
the central regions. 

for rotation and the velocity dispersion field presents a nearly circu- 
lar symmetric distribution (Paper III). We constructed axisymmet- 
ric /{EjL^) Jeans models as in the previous section, taking into 
account the apparent flattening of the galaxies in the outer regions, 
however, given the near-spherical shape in the central regions, the 
models are almost fully isotropic with / — It is sufficient to 

compare data and models along circular ringo The SAURON data 

^ The two-dimensional data versus model comparison for NGC 4552 was 
presented in Fig. Al of Paper IV. 



were averaged by computing a biweight mean and standard devi- 
ation jHoaglin. Mosteller. & Tukevlll983h of the values extracted 
from the Voronoi bins within radial ranges. The use of integral-field 
kinematics provides radial profiles of fi2 with very small scatter. 

If the galaxies were nearly isotropic in their central regions, 
as the Schwarzschild models show, the simple Jeans models, with 
only M/L as free parameter, should provide a good representation 
of the observed /i2 profiles. The comparison between the /i2 profiles 
of the data and the isotropic models (Fig. lA2t shows that indeed the 
models generally reproduce the data within the small uncertainties, 
and confirms the consistency of the observed galaxies with a nearly 
isotropic orbital distribution. The galaxies are of course not per- 
fectly isotropic: deviations between data and model do exist, partic- 
ularly in the outer regions of NGC 4486, and these small deviations 
from isotropy are consistent with the location of the galaxies below 
the isotropy line in the {V/a, e) diagram (Fig.|9]l. 



APPENDIX B: TESTING THE {V/a, e) DIAGRAM WITH 
ISOTROPIC MODELS 

For a quantitative use of the {V/a, e) diagram three questions have 
to be answered: (i) can one measure reliable values of the anisotropy 
using the (V/a, e) diagram from integral-field kinematics with lim- 
ited spatial coverage? (ii) Can one derive realistic anisotropics for 
models with an ellipticity that varies as a function of radius? (iii) 
what is the value for the parameter a (equation 1141 ) appropriate 
for realistic galaxy models? 

We tried to address some of these questions in Section |4!4l by 
comparing the anisotropy derived from the Schwarzschild models 
with the one measured with the {V/a, e) diagram. The good agree- 
ment between the two determinations showed that the answer to the 
first two questions is positive. These tests allowed us to verify the 
reliability of the {V/a, e) diagram for a wide range of anisotropics. 
However the approximate nature of the orbit-superposition models 
does not allow the test to be carried out very accurately. 

An alternative would be to use analytic models. The limitation 
of these tests is that they are easy to carry out only in the isotropic 
case. This was done by Binnev (2005) who also concluded with 
a positive answer to the first two questions above. In addition he 
measured a value of q < 0.2 for his test cases. However his adopted 
models were not entirely realistic, so the results cannot be necessar- 
ily generalised to real galaxies. 

In this section we revisit what was done by lBiIme^^ j2005f) . 
using the MGE models to approximate in detail the surface bright- 
ness of actual galaxies. For this test we selected five galaxies which 
represent some of the different types found in our SAURON sample 
and for which we derived a nearly edge-on inclination in Paper IV. 
The galaxy NGC 4486 is a giant elliptical with the largest ae in our 
sample. It shows no sign of rotation and has a nuclear cuspy core in 
the surface brightness profile. NGC 4552 and NGC 5831 are also 
classified as slow rotators and have central cuspy cores. NGC 3377 
and NGC 4621 are classified as fast rotators, are coreless and have 
disky isophotes. 

We computed the projected V and a fields for the five galax- 
ies under the assumption of axisymmetry and self-consistency (i.e. 
constant M/ L). In this case, when the distribution function is as- 
sumed to be a function f{E, Lz) of the two classical integrals of 
motion, the energy E and the angular momentum with respect to 
the symmetry axis, the non-centred second moment /i2 is uniquely 
defined by the density, once the inclination i and the M/L are cho- 
sen tLvnden-BelL, 196Z ; .Hunter,, 19 77). If in addition one assumes 
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isotropy of the velocity dispersion tensor, then the first moment fii 
of the velocity is also uniquely specified. Both moments /ii and fi2 
can be computed by solving the Jeans and Poisso n equations (e.g. 
ISatohlll98d : iBinnev. Davies . & Illine worthlll990h . In the case the 
density is described by the MGE parameterization the moment pi, 
projected on the sky plane, can be evaluated as a double integral us- 
ing equation (59) o f lEmsellem elaP ( Il994 . The projected moment 
jj,2 is evaluated as a single integral via equation (61) of the same 
paper (corrected for a typo as in Paper IV). We made the standard 
approximate identification V = fii and = ~ Mi- 

To be able to measure the parameter a the MGE models need 
to have constant ellipticity as a function of radius, as the adopted 
form of the {V/ a, e) diagram precisely applies only in that case. For 
this we took the MGE parameters as given in Paper IV but we set 
all the ellipticities e equal to a constant value e = 0.2, 0.4, . . . , 0.8. 
In this way the new MGE models have the same realistic major axis 
profile of the original galaxy, but have constant e. To evaluate the 
Jeans equations in a computationally efficient and accurate way up 
to large radii, and even for very flat models, we computed them on a 
grid linear in the logarithm of the elliptical radius and in the eccen- 
tric anomaly. This was done by defining a logarithmically-spaced 
radial grid Vj and then computing the moments at the coordinate 
positions {xj,k,yj.k) = [rj cos9k,rj{l — £)sin6'fe], for linearly 
spaced 9k values. The model was then re interpolated onto a fine 
Cartesian grid for the estimation of the V/a as done on the real 
integral-field data. 

In Fig. IB II we show the location of the five isotropic galaxy 
models on the (V/a, e) diagram. In the top plot the V/a was mea- 
sured from the model by integrating it within an ellipse of elliptic- 
ity e having mean radius i?* = a^/l — e = 3i?c, where a is the 
semi-major axis of the ellipse. For each different e the five models 
have similar V/a value and they lie along the theoretical isotropic 
line if one adopts oi ~ 0.20 ± 0.03. In the second panel the V/a 
was measured within an ellipse of mean radius J?* = Re, which 
closely resemble what we are doing on the real SAURON data. In 
this case the models lie slightly below the isotropic line, but the 
inferred anisotropy is still negligible 5 < 0.1. This shows that an el- 
lipse of mean radius Re is still sufficient for a reliable recovery of 
the V/a value of realistic galaxy models. This is not the case when 
the V/a is measured within an ellipse of mean radius i?* = Re/ 2, 
as in this case the values start being significantly underestimated 
{5 > 0.2, third panel). In the bottom panel the V/a was measured 
within circles of radius Re- The V/a values are now significantly 
below the isotropic line, particularly for very flattened models, and 
this shows that in the case of finite spatial coverage the V/a should 
be preferably evaluated along flattened ellipses than along circles. 

To test the reliability of the {V/a, e) diagram in a more real- 
istic case of variable ellipticity, we computed the (V/a, e) values 
inside l_Rc from isotropic models of the five test galaxies computed 
without changing the ellipticities of the Gaussians in the MGE mod- 
els of Paper IV. The ellipticity used is the ellipticity measured on the 
real galaxies. The result of this test is shown in the second panel of 
Fig. IB II It shows that the iV/a, e) values extracted from isotropic 
models of realistic galaxies, extracted within \Rc, do lie close to 
the theoretical isotropic line for oblate galaxies (when a ~ 0.2). 



APPENDIX C: MONTE CARLO SIMULATION 

In Section [5^ we showed that the location occupied by the fast 
rotating galaxies in the SAURON representative sample is consis- 
tent with the one expected according to equation l |19i l. It is however 
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Figure Bl. V/cr of isotropic galaxy models. Top plot: the V/cr was mea- 
sured from the model by integrating it within an ellipse of ellipticity e and 
mean radius Ri, = a^/l — e = 3Rc, where a is the semi-major axis of 
the ellipse. Second panel: the V/ u was measured within an ellipse of mean 
radius Ri, = R^. The NGC numbers refer to the {V/u, e) values measured 
on the isotropic models with variable ellipticity. Third panel: V/cr within an 
ellipse of mean radius R^, = Re/2. Bottom panel: the V/a was measured 
within circles of radius Rg. This plot shows that V/cr is better recovered 
along ellipses than along circles of the same area. Moreover the V/ cr values 
appear to be well recovered up to an ellipse of mean radius Rc . 
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Figure CI. Logarithmically spaced contours of the likelihood L{ms, crs) 
as a function of the assumed slope and spread of the 5 — e coiTelation. The 
maximum likelihood values are mg =0.61 and ag = 0.11. 




Figure C2. Same as in Fig. |9] with overlaid with the magenta lines the 
contour of the probability P' [(V'/a-)obsi ™(5i ""is]. for the best fitting pa- 
rameters of Fig. ICll The contours are meant to describe the distribution of 
the fast rotating galaxies (blue labels). 

unlikely that galaxies will follow a correlation with zero scatter. 
Moreover the previous discussion does not make use of the infor- 
mation contained in the whole distribution of observed galaxies in 
the (V/a, e) diagram. As we do not know the inclination for all the 
galaxies in the sample, the only way to verify if the observed mea- 
surements imply a correlation between shape and anisotropy is to 
resort to Monte Carlo simulations. 

We saw in Section 15.11 that the fast rotators are consistent 
with being a population of nearly oblate galaxies. In principle, un- 




E 



Figure C3. Monte Carlo simulations of the [V/a, e) diagram. Top panel: 
the galaxies follow a relation S = 0.6 eintr with added Gaussian intrinsic 
scatter. Second Panel: same as in top panel, but galaxies with <5 > 0.8 eintr 
are removed. Third Panel: at every Eintr galaxies are uniformly distributed 
in the interval 5 = [0, O.Seintr]. Bottom Panel: galaxies are on average 
isotropic 5 = 0, with some added Gaussian intrinsic scatter. Only the second 
panel resembles the observed distribution of the fast-rotators. 
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der the oblate hypothesis, it is possible to invert the observed el- 
lipticity distribution to derive the intrinsic shape distribution, un- 
der the assumption or random orientations. However the SAURON 
sample was explicitly selected to cover a nearly uniform range 
of ellipticities (Paper II), and this means that: (i) one cannot as- 
sume random orientations and recover the intrinsic distribution, but 
also that (ii) the observed distribution on the {V/a,e) diagram 
will be weakly dependent on the assumed ellipticity distribution, 
once the same selection criteria are applied. For this reason, while 
still assuming the fast rotators are oblate, we will adopt the para- 
metric form for the distribution of in trinsic ellipticities derived by 
iLambas. Maddox. & Lovedavl ( Il992[) for a sample of ~ 5, 000 SO 
galaxies: 



L{ms,as) = J| P' l{V/a) 



(C5) 



'(Sintr 




- 'n-.T'l 'f0<eint.<0.65 ^^^^ 
otherwise, 



with eo = 0.41 and = 0.24, and where we set an upper limit on 
Eintr based on the maximum ellipticity observed in our sample. 

For any given ellipticity we assume a simple linear relation 
with slope ms to determine the anisotropy, as found in Section l46l 
but we allow for an intrinsic (Gaussian) scatter around the relation. 
The distribution of the anisotropy 5 at any given ellipticity is then 
assumed to be: 



exp 







otherwise, 



where the upper limit for S comes from the condition (V/a)'^ > 
in equation il4l and {ms,f7s) are two free parameters. 

Given Sintr and 5 we use the inverse of equation l ll4t to 
compute the edge-on V/a value. Finally with equations il2ll3t 
we project the {V/a, eintr) pair of values to get the observed 
[(y/cr)obs, e] and {V/a)obs values, assuming a uniform distribu- 
tion on the sphere of the viewing angles (i. e. a uniform probabil- 
ity distribution in sin t). This process, repeated for a large number 
~ 10** of random realisations of the of the eintr variable, is used to 
compute an histogram. This constitute a numerical estimate of the 
probability P[{V/a)ohs,£', ms, ag] of observing a galaxy within a 
certain range of values on the {V/a, e) diagram, for any given set 
of our two model parameter {ms, as). 

One complication now arises because of the SAURON selection 
criteria. We need to consider the fact that the value of Eobsj was 
specifically chosen for each sample galaxy, and does not come from 
a random selection. For this reason we define a new probability 



P'[{V/a)ohB,£;ms,as] = 

P[{V/a)ohs,e;ms,as] 

/o°° P[{V/a)ohs,e;ms,as] d{V/a)ob 

in such a way that 



P'[{V/a)ohs,e;ms,as] d{V/a)obs = 1, 



(C3) 



(C4) 



or in other words the marginalised probability of observing a galaxy 
with any given e is constant. What varies is only the probability of 
observing a certain {V/a)ohs value at the chosen e. 

To find the best fitting parameters of our simple model we then 
define a likelihood function, which gives the probability (neglecting 
the effect of measurements errors) of observing the given set of M 
independent measurements [(V/(T)obs,m, eobs,m], for the assumed 
model parameters {'m.s,as): 



We determine the likelihood only for the fast rotators of the 
SAURON sample of 48 galaxies, and excluding the galaxies 
NGC 2685, NGC 3156, as they lie in regions of low-number statis- 
tics and may not follow the general distribution of the other galax- 
ies, and NGC 4473, for the reasons explained in Section [5^ This 
leads to a sample of M — 36 galaxies. The contours of the like- 
lihood function l |C5t are presented in Fig. ICll and show that the 
maximum likelihood parameters are ms = 0.6 and as — 0.1. 
Given the very different approach adopted for their determination, 
the agreement with equation l |19t provides a strong confirmation of 
the observed trend of anisotropy versus ellipticity. The contours of 
the probability P' for the set of best fitting parameters are shown 
in Fig. |C2| Although the contours follow the general distribution of 
the fast-rotators, they also display a tail with non-negligible prob- 
ability at low values of V/a, where we do not observe any galaxy. 
This effect is also visualized in the top panel of Fig. |C3l where we 
plot the distribution of values on the {V/a, e) diagram, for a simu- 
lated set of 300 galaxies. A number of galaxies are predicted by this 
model to fall below the magenta line. 

After some experimentation we found that the simplest way 
to reproduce the observations consists of enforcing the upper limit 
on the anisotropy 5 < 0.8 eintr in equation l |C2K The result of this 
simulation is shown in the second panel of Fig. IC3I and qualita- 
tively resembles the appearance of the observed SAURON {V/a, e) 
diagram. The distribution of the observed galaxies in the {V/a, e) 
diagram suggests that the above condition constitutes a physical 
zone of avoidance for the existence of the fast-rotators. This zone 
of avoidance may be related to the stability of axisymmetric stel- 
lar systems with varying degrees of anisotropy and would deserve 
further investigations. 

In Fig. |C3| we present two additional Monte Carlo simulation, 
under simple assumptions for the galaxy anisotropy. In the third 
panel we assume, at every eintr, an uniform distribution for galaxies 
in the interval 5 — [0, 0.8 eintr]- This distribution is unable to repro- 
duce the observed distribution of the fast-rotators on the {V/a, e) 
diagram. Finally, in the bottom panel we assume that galaxies are on 
average isotropic <5 = 0, but we allow for the same Gaussian intrin- 
sic scatter as = 0.1 as in the top panel. This distribution is strongly 
excluded by the observations. A comparison with the {V/a, e) dis- 
tribution for a much larger and unbiased sample of galaxies would 
be required to further confirm these results. 
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